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ABSTRACT 


A  procedure  co  approximate  data  given  at  arbitrary  intervals  in 
two-,  three-  and  four  dimensions  using  polynominal  expressions  Is 
described.  The  programing  of  the  problem  is  explained  in  each  case 
and  a  user  manual  is  given  for  softvare  implemented  on  the  TPL 
Hewlett-Packard  computer  system.  The  method,  which  is  general,  was 
derived  and  has  been  applied  to  represent  Che  calibration  of  flow 
probes  which  have  multiple  sensors. 


1.  INTRODUCTION 


This  report  describes  algorithms  to  approximate  data  patterns  such 
as  those  which  occur,  for  Instance,  when  pneumatic-velocity  probes  are 
calibrated.  The  functional  value  (Y)  can  depend  on  either  one  (XI), 
two  (XI  and  X2)  or  even  three  (XI,  X2  and  X3)  parameters.  The  approxi¬ 
mation  Y  -  f (XI)  will  be  referred  to  as  the  two  dimensional,  the  approxi¬ 
mation  Y  -  f (XI,  X2)  as  the  three  dimensional  and  the  approximation 
Y  «  f(Xl,  X2,  X3)  as  the  four  dimensional  approximation.  These  three 
options  meet  the  requirements  for  probes  used  in  the  Turbopropulsion 
Laboratory. 

Using  the  least-squares  criterion  to  obtain  the  coefficients  in 
assumed  polynomlnals  leads  to  a  system  of  linear  equations,  which  are,  in 
principle  easily  solved.  Numerical  problems  may  arise  however,  since 
the  software,  as  described  herein  and  as  Implemented  in  the  UPS 
Turbopropulsion  Laboratory  Hewlett  Packard  21MX  computer,  uses  32-bit 
real  constants. 

The  author  wishes  to  express  his  thanks  to  Professor  Ray  Shreeve, 
whose  constant  and  critical  Interest  was  a  crucial  help  to  solve  this 
mathematical  problem  and  link  it  to  engineering  application. 
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2.  TWO  DIMENSIONAL  APPROXIMATION 


2.1.  Problem: 


A  data  sec  of  NPNTS1  data  points  Is  given,  where  7  depends  on 
parameter  XI.  The  data  pattern  Is  to  be  approximated  by  a  function 
7  *  f (XI) ,  so  chat  the  error  between  data  points  and  analytically 
determined  points  Is  lowest. 


2.2.  Approach: 


Polynomlnals  ara  commonly  uaad  In  ordar  co  approximate  data 
patterns.  T  »  f (XI)  is  a  function  which  approximates  the  data  value, 
V,  at  each  value  of  the  variable  XI,  we  look  for  our  expression  of 
the  form 

r  -  c  ♦  c  *xi  ♦  c  *xi2  ♦  ...  ♦  c  *xi<L”1> 

12  3  L 


or 


Y 


-  fc  C  *Xi 
i*i  i 


<  i-i) 


(2.1.) 


in  which  tha  coefficients  are  to  be  determined  by  the  method  of 
least  squares.  As  Ref .  1  shows  the  least  squares  criterion  leads  to 
a  linear  equation  system  that  can  easily  be  solved.  We  define  the  error 


MPNTS1  2 

*  •  E  tf(Xi  >-Y  1 
r»l  r  r 

(2.2.) 

where  the  index  r  denotes  the  individual  data  point.  Using  aquation 
(2.1.)  R  becomes 


*  - 


HP  NT 


r«i  i 


♦  C  *X1  ♦  C  «Xi  ♦ 

2  r  2  r 


.  C  .Xl"-1'  -  Y  ,2 
L  r  r 
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2.3.  Solution: 


R  depends  on  Che  selection  of  the  coefficients  C. ,...,C  . 

X  L 

In  order  to  minimize  the  erTor,  R  is  differentiated  with  respect 
to  C^,...,CL  and  the  partial  derivatives  are  set  to  zero.  Thus 


3C 

x 


or 


3  R  NPNTSl  2  (L-l) 

- —  *  £  2*tC  +C  *Xi  +C  SXi  +  .  .  .  ♦  C  XI  -Y  ]* 

3C  r»  1  1  2  r  3  r  lr  r 

i 


^  IC  *C  SXI  +C  SXI*"*  ...  +  C  Xi<L_l>-Y  1 
12r3r  L  r  r 


Assuming  that  the  sunmatlons  extend  over  all  the  data  points,  Z 

NPNTSl 

should  be  understood  to  mean  E  .  Performing  the  different 

r«i 

tlations  and  rearranging  the  equations,  we  get 


NPNTS1SC  +£xi  sc  ♦  £xiasc  ♦  ...  ♦  £xi<l"usc  -  £y 
x  r  2  r  3  L  r 

Hxi  *c  *£xi2sc  ♦  £xi\c  ♦  ...  ♦ExiU  sc  -Ey  sxi 

r  i  r  2  r  3  r  I  rr 

sc  ♦  £Ixi  sc  +27x1  sc  ♦  ...  *£xi  L  l>sc  *Cy  sxi2 

r*  rg  r  3  r  Lrr 


27xi 'L'1 'sc  ♦Z?xiLsc  ♦  Exi<L+nsc  ♦.  .♦£xi<2L”2>sc  -£y  ,xil“ 

r  *  re  r  3  r  Lrr 
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In  matrix  notation 


NPNTSi 


EXI 

r 

2 

EXi 


L-i 

EXi 
'  r 


2 

L-l 

Exi 

Exi 

...  Exi 

C 

Ey 

r 

r 

r 

i 

r 

2 

3 

L 

EXi 

Exi 

.  . .  EXi 

c 

EY  *Xi 

r 

r 

r 

2 

r  r 

3 

4 

L+l 

* 

9 

2 

nxi 

EXi 

...  EXi 

C 

EY  *Xi 

r 

r 

r 

3 

r  r 

L 

L  +  i 

2L-2 

L-l 

EXi 

EXI 

...  EXI 

C 

EY  *X1 
r  r  f 

r 

r 

r  ' 

^  L  ' 

* 

9 

X*  (2.3.) 

f5*  ...  Systen  Matrix 

...  Coefficients  vector 
E«  ...  Right  hand  side  vector 


8 


2.4.  Equation  System 
2.4.1.  System  Matrix  A 


Since  all  Matrix  elements  on  diagonal  lines  running  up  from 
left  to  right  (2.  order,  as  defined  in  Appendix  A)  are  identical,  the 
elements  are  renamed  for  simplification  as  follows: 

<2,  *  *f 

ar-  V  a* 

aJ  * 


ai  * 

*  a 

m  -• 

'  *2,"  "  V 

'  aL-t,J  *  — 

*  «JE,AW  *  *2.L 

?L+2 

d2l-l 

X(*  LjJ,  = 

*  ait  4W  *  Gs,t. 

Q 

so  that  the  set  of  elements,  a^,  can  be  written  as 
tmmt  . 

ak  »  27*7  /  ... 2L-! 

r*»  1  '  (2.4. 

This  operation  is  programmed  in  REAL  FUNCTION  S2  (NPNTS1, IP0WR1, 17 ) . 
The  data  XI  and  T  are  known  to  this  function  through  a  cotmon  block 
named  DTA2.  So  only  NPNTSI,  EP0WR1  (-  k-1)  and  IT  (-0)  have  to  be 
passed  to  the  function,  and  the  value  of  a^  is  returned  through  the 
function  name  S2. 
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SUBROUTINE  MAT 2  presets  the  system  metrlx  in  the  following  way: 


<*u 

a,t 

av 

<%I 

°xi 

<** 

•••  •••  .  A»** 

S  c/  c/ 

0  0  • 

*'7 

Qv  ••• 

°.L 

<*u 

Qxx  ••• 

<*ZL 

Oj, 

Off  •  • » 

sea 

•  •  • 

au 

#*•  • • • 

••• 

<*a. 

2.4.2.  Right  hand  aids  vector  B 


Th«  elements  of  the  right  hand  slda  vector  B  can  be  written  as 

NTKIJ) 

m-xi/;  k.  ... jL 


To  calculate  b^,  again  REAL  FUNCTION  S2  la  used.  Data  XI  and  Y  are 
available  through  COMMON  block  DTA2.  NPNTS1,  IPOWR1  (-k-1)  and  IT  (-1) 


have  to  be  passed  to  the  function,  and  the  value  of  b^  is  returned 
through  the  function  name  S2. 
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2.5.  Software: 


The  software  to  compute  the  coefficients  for  a  two  dimensional 
approximation  is  described  in  Appendix  B  and  is  implemented  in  the 
TPL  HP-21MX  computer  system 

To  work  correctly  with  these  program  modules,  the  user  must 

conform  to  the  following  conventions: 

i)  Provide  the  data  in  two  arrays  (Type:  REAL)  of  256 
elements  through  a  COMMON  block,  named  DTA2. 

COMMON  /  DTA2  /  XI, 7 
REAL  XI (256) ,7(256) 

11)  Dimension  an  array  (Type:  REAL)  of  7  (seven)  elements 
to  contain  the  coefficients. 

REAL  COE? (7) 

Hi)  Define  the  parameters  NPNTS1,  L  and  EPRINT  (Type:  all 
INTEGER) 

NPNTS1  ...  #  of  data  points 

L  <  NPNTS1  <  256 

L  ...  (desired  order  of  polynomlnal)  +  1 

1  <  L  <  7 

IP HINT  . . .  controls  quantity  of  print  out 

2  ...  Print  system  matrix  and  right 

hand  side  vector  before  and 
after  Gauss  Jordan  Elimination 
1  ...  Print  equation  system  after 

Gauss  Jordan  Elimination 
<0  ...  No  print  out 

>0  ...  Print  equation  (1,1)  with  the 

actual  paramatars 

iv)  Whan  loading  a  program,  that  uses  the  subroutine  MAT2, 
the  binary  library  file  has  to  be  searched  for  externals. 

Suppose,  the  source  file  &USER: :26  contains  a  user 
program  named  USER.  This  program  calls  MAT2.  After 
the  compilation,  tha  relocatable  binary  file  for  this 
user  program  is  XUSER: : 26.  To  load  tha  program,  the 
following  procedure  is  recommended: 


:RUf LOADR 


and  che  loader  program  will  respond 
/LOADR:  RE.ZUSER: :26 

where  che  underlined  informacion  already  is  Che  user's 
Input.  This  causes  Cha  loader  Co  load  all  program 
nodules  of  che  user  program;  a  load  map  is  listed  on 
che  eernlnal.  Upon  completion  che  loader  prompts 

/LOADR:  MS.ZTPLBL: : 26 

where  che  underlined  information  already  is  che  user's 
Input.  Now  che  loader  conducts  a  search  for  all  unsatis¬ 
fied  externals  of  che  user  program  in  the  binary  lib¬ 
rary  file  ZTPLBL: :26.  Since  some  library  programs 
have  externals  themselves  che  search  has  Co  be  repeated 
(MS  . . .  multiple  search)  as  many  time  as  is  necessary 
for  all  externals  Co  be  satisfied.  The  loader  prompts 

/LOADR:  END 

where  che  underlined  information  is  the  user's  input. 

The  loader  now  loads  all  system  programs,  outputs  a 
load  map  and  generates  the  program.  Upon  completion, 
che  loader  outputs  a  ready  message 

/LOADR: USER  READY  AT  1:28  PM  WED.,  10  SEPT,  1980 
/LOADR:  END 

Now  che  program  can  be  run 
:RP.PSER 

Again,  che  underlined  informacion  is  che  user's  input. 

If  all  these  requirements  are  met,  the  correct  call  for  che  sub¬ 
routine  is: 


CALL  MAT2  (NPNTS1 ,L,C0EF, IPRINT) 

Upon  completion,  che  array  COEP  contains  the  coefficients  (COEPCD-C^ 
, . . . ,C0EF(L)*C^) .  Externals  used  by  MAT2  are:  AB2,  DTA2,  S2  under 
no  circumstances  may  the  user  use  any  of  these  nsmes  for  modules  of 
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his  own  user  progrsa.  In  some  cssss  ths  program  may  not  be  able  to 
perform  a  Gauss  Jordan  Elimination  to  the  system  matrix  and  the  right 
hand  sida  vector  and  thus  causa  the  program  to  stop.  If  this  happens, 
an  erre*  message  is  displayed. 

It  is  highly  encouraged,  to  oae  the  system  function  FNP  to 
calculate  the  value  of  a  polynominal  at  one  specified  XI  rather  than  a 
loop  such  as  the  following: 

S-.O 

DO  08  I-l.L.l 
08  S«S+X**(I-1)*C0EP(I) 

YCALOS 

To  call  FNP  is  easier  and  faster,  because  FHP  is  a  prograamad  Horner 
Schema.  The  recommended  call  is  than; 


Ll-L-1 

YCALOFHP(COEF,X,U) 


2.6.  Sample  User  Program 


2.6.1.  Listing 

PACE  0001  PTH,  12  >  40  PH  THU.,  18  SEP.,  1980 


0001 

0002 

tm 

ooos 

0006 

0007 

0008 

0009 

0010 

uu 

0013 

0014 

001S 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

124 

.025 

0026 

S3I3 

0029 

0030 

Ji3i 

0033 

0034 

Soil 

0037 

0038 

0039 

9040 

0041 

•  042 
0043 
0044 
0045 
S  @44 
0047 

•  048 
“  *149 

JS® 

3311 

0053 

0054 

im 

0057 

0058 

Mil 

0061 

0062 

ml 

0065 

0*66 

ms 

0069 


FTN4,L 


PROGRAM  DEM02  <3, 99) 


THIS  IS  A  DEMONSTAT1QN  PROGRAM  AND  IT  SH0W8  THE  CORRECT  USE 
OF  THE  TPL  BINARY  LIBRARY  (ON  TYPE  6  FILE  XTPLBL)  USING 

smi§LT  ft°hI  th£ 


i,  DEMONSTRATE  THE  USE  OF  THE  TPL  BINARY  LIBRARY. 
COMMON  /  AFLD  /  A 


REAL  A<256> 


NTECER  NOLF , NOCR  <  2  > , 1 CLR  ( 3  > 

DATA  NOLF  /006537B/ 

DATA  NOCR  /0©iQ33B, 040433B/ 
DATA  I CLR  /015524B?0iSSl5B,! 
-  /3. 14159$/ 


DATA  PI 


006S37B/ 


c  1.1  fo&SW/J-3  *Rfi  IB?  IsSSiS&'W  AND  ME ’LL  SEE  H  „„ 

tOU  TO  USE  THE  MARVELLOUS  TPL  BINARY"A2/"  LIBRARY!  YOUR  INTEREST 
•WILL  BE  GREATLY  AWARDED  BY  EASIER  PROGRAMMING. ■///> 

US  mm  '.a.fFIH’ar  wmwiin<J»W(«>-t 

m  foSSSt  '<• 

m  V-  .EWrfHqiWUG.  AREAS  AND  DRAW  AXES") 

108  FORMAT  <9X""6X"XMIN  "6X"XNAX  "6X"YHIN 
*"6X"YMAX"/9X.4<""10X"  ■)/ 

*-  HP9872",4<"  "F10, 3" 


*"  QSEKZ 

109  FORMAT^1 

110  FORMAT  ( 

Hi  ESIffiST  I 

*•) 

113  FORMAT  ( 


V  uim  tlnW  L AsuEF I  NED  AND  AXES  DRAWN!"/) 

DRAWING  DATA  POINTS  INTO  COORDINATES  SYSTEM") 

HalculaHnS  ARCURV!E/FIT  THROUGH  THE  DATA  POINTS 


113  FORMAT  (•  ENTER  DEGREE  OF  POLYNOMIAL  (NORDER)  TO  FIT  THRO 
U4*$5rNATE<  "01EN?ER  fpSlNT  "2A2> 

«f  FORMAT  V' 

!»  am  v- 

C  149  F0RMAt|W>DEM02  stop 


GET  THE  LU  OF  THE  TERMINAL. 


iJlTE^ft^lOi)  NOLF 


GENERATE  DATA  POINTS  SCATTERED  AROUND  A  POLYNOMIAL. 


0073 


call.  Mtfc’  my 

MITE. (LI ,  102) 
(START  -  -6. 
(STOP  •  *2. 
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PAGE  0002  DEH02  12>40  PH  THU.,  18  SEP.,  1980 


0078 

0079 


0083 

0084 

008S 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

009S 

0096 

0097 

"*»98 

J99 

0100 

0101 

0102 

0103 

0104 

0105 

tlti 

0108 

0109 

0110 

0111 

".ill 

0116 

0117 

0118 

0119 

0120 


ii 
0 

o.» 

0128 

0129 

Q4  30 

0 
0 
0 

0l3S 

•oil* 

Sil? 

0140 

0141 

na 

0144 

0145 

im 


OX  «  (XST0P-XSTART)/(NPNTS1-1) 

COEF ( 2 )  -  -0.500000 
CQEF<3>  -  -0.000340 


CuEF < 4 )  »  ♦0.300000 
COEF(S)  -  ♦O.OSSOOfS 
A  -  <2. OftPDOO . 03125 

5?<?i 

Y< 1 )  »  FNP(C0EF,X1(X) ,NORDER)  ♦ 

01  WRITE  (LI,  103)  I, XI (I) ,Y(I) ,  NOLF 


0.2S*SXN(X1(X)/A) 


DUHHY  •  T<1+1> 


Y(X4l)  -  Y 
02  Y(I)  -  D 

Miff  St!; 


■  Y<  I > 

-  DUHHY 
LI,  149) 
LI,  104) 


< ICLR , X"1 ,3,1) 


INITIALIZE  THE  PLOTTER. 


WRITE  (LI,  10S) 

OH  ICLR 
WRITE  (LI,  106) 


DEFINE  PLOTTER  AND  USER  AREAS)  DRAW  AXES. 


WRlti  (LI.  10*) 

XPHIN  -  2. 

XPHAX  -  12 . 

YPHIN  -  12. 

YPHAX  ■  22 . 

ass  : 

YUHIN  ■  -5 . 

YUHAX  >10. 

WRJTEx(LI,.lJJ>  XPHIH, XPHAX, YPHIN, YPHAX, XUHIN,XUHAX, YUHIN, YUHAX 

XA  -  (XPHAX  -XPHIN  >/(XUHAX-XUHXN) 

XB  -  ( XPHIN*XUHAX-XPHAX*XUHIN ) / ( XUHAX-XUHIN) 

a>H*y :  WW  -X'HIN 

YA  -  (YPHAX  -YPHIN  >/( YUHAX- YUHIN) 

YB  »  ( YPHINBYUHAX-YPHAXBYUHIN  > / ( YUHAX- YUHIN  > 

YL  ■  (YPHAX  -YPHIN  > 

CALL  8ETSH  (113,1.) 

KttKil 
Sift  St!;  151! 


PLOT  DATA  POINTS. 


•  liiil* . 

WRITE  (LI,  110) 


PAGE  0003  DEM02  12»40  PH  THU.,  18  SEP.,  1980 


XPLOT  -  XI < I ) 4XA+XB 
YPLOT  ■  Y(I)*YA+YB 
_  CALL  PLOT  <XPLOT, YPLOT, 2> 
03  CALL  SYHBL  (1) 

WRITE  (LX,  149)  XCLR 
WRITE  (LI,  ill) 


CALCULATE  CURVE  FIT  THROUGH  DATA  POINTS. 


"M  It I: 


04  UR 
READ 
URITf 


i>2) 

113)  NOCR 
»)  NORDER 
149)  ICLR 
.LT.  0 


URITt  1|4)  NOCR 

(ICLR7l-i,2,l> 


OR.  NORDER  .GT.  6  )  GO  TO  04 


IF  (XPRlNt.GE.O)  WRITE  (LI,  11S> 
IF  (IVRINT.LT.O)  WRITE  (LI,  116) 


SiSi 


PLOT  POLYNOMIAL  FIT. 


WRITE  (LI,  117) 

CALL  SETS*  (113.2.) 

DO  89  I-1.NPNTS1.1 
XPLOT  -  Xl(X)tXA+XB 

ifl?I  .  e5 . 57  JffiFrfW1  \  dmw X1 

09  IF  (I.GT.l)  CALL  PLOT  (XPLOT, YPLOT, 3) 
CALL  SET8M  (US, DUMMY) 

tlifi  ■<«-,  mi ia ■* 


TERMINATE  GRAPHICS. 


CALL  STOPG 
|TJP  7777 


FTN4  COMPILER  t  HP92060-16092  REV.  1926  (790430) 

**  NO  WARNINGS  **  NO  ERRORS  **  PROGRAM  •  01409  COMMON  «  00000 
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2.6.2.  Load  aa 


0EHQ2  10042  12642 
FNP  12643  12746 


2.6.3.  Results 


Affeet  of  IP  RIOT  on  quantity  oi 


1)  printed  output 

System  Matrix  A  and  vector  B  before  Gauss  Jordan  Elimation 


1  .5lQE+0^ 

2  - . 102E+03 

3  . 487E+03 

4  -.211E+04 

5  .  104E+0S 

i 


102E+03 

407E+O3 

211E+04 

ttn:n 


487E+0’ 


,  21 1E+04 


21 1E+04  .104E+05 

i 04E+05  - . 525E+05 

s$ii:§i  --.ntm 

3  4 


104E+oI  -.S78E+02 

1111:81  -:II!i:8i 
mM  -.IHBI1 

s 


Equation  system  after  Gauss  Jordan  Elimination 


i  . lOOE+Oi  .000E+0 


OOOE+OO  . OOOE+OO  . OOE+Ol 


970E+00 


2  . OOOE+OO  .100E+01  .OOOE+OO  .000E+00  . OOOE+OO  -.416E+00 

3  .OOOE+OO  .OOOE+OO  .100E+01  .OOOE+OO  .OOOE+OO  -.1A2E-01 

4  .OOOE+OO  .OOOE+OO  .OOOE+OO  . lOOE+Ol  .OOOE+OO  .27SE+qO 

5  .OOOE+OO  .OOOE+OO  . OOOE+OO  . OOOE+OO  . iOOE+Oi  .514E-01 

1  2  3  4  5 


IPWNT-2 


Coefficients  COEF(I) 

I  -  1  - . 970E+00 
1-2  -.416E+00 
1*3  - . 162E-01 
1-4  . 275E+00 
I  «  5  .S14E-01 


Equation  system  after  Gauss  Jordan  Elimination 


1  . 100E+01  .OOOE+OO 

2  .OOOE+OO  . 1Q0E+01 

3  .OOOE+OO  .OOOE+OO 

4  . @00E+0Q  .OOOE+OO 

5  .OOOE+OO  .OOOE+OO 

1  2 


Coefficients  COEF(I> 

\ :  i  -Kmn 

i :  3 

I  -  5  . 514E—01 


OOOE+ol  .OOOE+OO  . OOOE+ol  -.970E+00 
OOOE+OO  .OOOE+OO  .OOOE+OO  -.41AE+00 
100E+01  .OOOE+OO  .OOOE+OO  - . 162E-01 

8!!i:°o8  -mtn  lissits;  .iniTi 

3  4  5 


I PRINT-  1 


Coefficients  COEF(Z) 
I  -  1  - . 970E+00 

\  :  l  -.tlSi-Sl 

in  . ini-Si 


«wr<  I 
>2 
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3.  THREE  DIMENSIONAL  APPROXIMATION 


3.1.  Problem: 


A  data  sat  of  NPNTS1*NPNTS2  data  points  is  given,  where  Y 
depends  on  the  parameters  XI  and  X2.  The  data  pattern  is  to  be 
approximated  by  a  function  Y»f  (XI,  X2) ,  so  that  the  error  bet*  een 
data  points  and  analytically  determined  points  is  lowest. 


3.2.  Approach: 


As  for  Cha  two  dimensional  approximation,  a  poly-nominal,  now 
with  two  independent  variables  XI  and  X2,  is  used. 

y-  fcxi.  *2) 

Y-  ca-xj£+  ...  *  C,M-X2 + 

*  ...  *  C2MXz""]-XI  * 

*  CCgi  "*  CfgGL  .  £33^2.  t  . , ,  *  Cjfj  X2  3  'Xt  * 


/  (fif-  tj 

t  C<3t,i  CuXZ-> CUXZ  ■>  ...  *0^X2  3-X! 


L  M 

E 

c*  / 


Y-T,{  Z.C,iX2(i">]-Xl‘"> 
y- '  / 


(3.1.) 


The  data  arrangement  for  the  three  dimensional  approximation 
already  yields  to  the  application  of  this  method:  at  each  constant 
Xl  (i.e.:  r  ■  constant,  too)  there  is  a  data  set  for  various  X2 . 

This  is  the  case  for  example,  when  a  pneumatic-velocity  probe  is 
calibrated.  At  each  Mach  number,  the  probe  is  balanced  in  yaw  and  the 
pitch  angle  is  varied  over  Che  range  of  application.  The  probe  out¬ 
put  (each  of  two  separate  non-dimensional  pressure  differences)  is  a 
function  of  the  Mach  number  and  the  pitch  angle. 
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The  least  squares  criterion  becomes 


R- 


lews/  -  HfMnt  —2  1 

Z-f  r  Cfrn„'XZHi-  Yn  J  ] 


0.2.) 


where  the  Indices  r  and  s  denote  the  Individual  data  points.  Using 
equation  (2.1.)  R  becomes 


r.  z-f  n  c  c.  ■>  .. 

r-y 

+  (CZ,*  — 

*£Cj,4CjiXZft+  ... 


.. .  cm>C  - 

■  *  c*h  *C ')■ 


The  term  In  the  [  J  -  bracket  shall  be  called  B. 


3.3.  Solution: 


To  minimize  the  error,  R  is  partially  differentiated  to  the 
coefficients  and  then  the  partial  derivatives  are  set  to  zero 


a*  i 


/  *  /  ..  4.  » 


ns-M  i  ■  o 

dCj  r-t  f-t 


Assuming,  that  the  summations  extend  over  all  data  points,  El 


should  be  understood  to  mean  j  ,«  jT  T7  Performing  the  differen- 

/  t 

tlatlons  and  rearranging  the  equations  in  matrix  notation  (which, 
while  lengthy,  is  logically  straight  forward),  we  get,  in  matrix 


notation 


EE  x/K 

ZX/ljXZ, 

EE  xi:xi 

EEx/vtZ' 

,  _ .  u-> 

ZZx^xz 

TZuxz, 

LLxin 

2  £*<,•«„  - 

M>vC 

EEjw* 

EE*/’ ■*£ 

ee*u: 

T.T.XI  m 

a  rr 

EE*'*-*2* 

! 

EEjt  /„•«* 

EE*£  *z 

ZZxS.az* 

... 

EE« 

EE«-«r 

rx  rr 

EE  */*«,"" 

ZX<-»* 

... 

EE^rC 

2*M 

S - 

(L-!)*  fi*l 

~  1 

EE  x&C 

ci>C 

£EV*„ 

... 

zzx£-*C ' 

L£*i-*z2 

< 

... 

ez'£«. 

22*CV. 

••  •  - 

S2« 

ZZuxfJ 

22<*«' 

t#« 

22<-*C 

EL^v«r 

2  &<C 

22<-«vl 

zjx-< 

»..  rsx<! 

#»» 

r2X*2„ 

ZZx£-*zl 

•  #  # 

se^*>2; 

2sx<1 
— - 1 

i 

' 

1  -1 

*»• 

ss>C«C 

Z2K< 

»•# 

_ .  »»•« 

LE 

Z7>I>r  ZLK'xz 


3.4.  Equation  System: 
3.4.1.  System  Matrix  A 


This  system  matrix  needs  some  examination.  A  crucial  step 


is  the  introduction  oi  submatrices  1.  class  into  the  system  matrix. 


4 

■4  ,L 

A  - 

A  I, 

2. 

The  asterisk  *  denotes  a 

submatrix  1.  class  where 

t+NVSHNPNTSZ 

EE*2„ 

ee 

<• 

LHxza 

EI< 

ET-C 

ILxir 

XZjcC 

'  Z£<"4 

HY>xin 

ZZxi„n“ 

<* 

DI 

EE  w*< 

EE^UV 

ZZxuC 

£2vC 

z&yC 

££*/» 

ZZxt*xza 

EE 

Z£*'X 

EE« 

EE 

EEJ$C 

ZHx£*2? 

Ex<«rJ 

etc. 
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Now  ic  can  seem  Chat  all  submatrices  on  the  same  diagonal  bands 


2.  order  are  Identical.  The  submatrices  A  subsequently  are  renamed 

/*  -  * 

/  *  A  // 


A*  * 

a* 

■  K 

mU 

A3  9 

■A* 

- 

A\  ' 

f 

K 

*  /A 

Ac  * 

K, 

■  <z 

-/t 

A* 

L+> 

A* 

*** 

/J 

« 

iT 


J  ~ 


Z  H 

L-f 


l 

'4,  4 


A 


M-l 


A 


can  be  written 

as 

xz*c>< 

EE<V' 

nr»  W  / 

Elxi„  xfc,. 

ss«r«; 

*■*  «  . 

f* 

'  K 


I 

s  I 

Like  the  system  matrix  for  the  two  dimensional  approximation  this 
submatrlx  1.  class  is  not  only  symmetrical  to  its  main  diagonal  line 
1.  order,  but  also  the  elements  on  each  diagonal  line  2.  order  are 
identical.  The  elements  a^,^  of  Che  submatrix  A^  therefore  are 


■■  ■  *  -  -  <»- 


V1'1 
'•NT  •• 
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;/  *  ak;i/ 

a*)L  -  aK;a  *  a«;tz 


aKf 3  *  a  A;  7/ 

"  *  a«}/3 

• 

ak\  H  *  aKf  Ml 

s  2  *  a«',H-2  3  *  *“ 

a*zin-i 

®  tn 

'  i  * 

‘  &K;  3  h-t 

-  aH . hi  *  ... 

*  O*  •  f  H-4 

-  <*#;  or  h 

G'k'ZH-/ 


Equation  (3.4.)  therefore  is  the  fundamental  step  towards  the  pro¬ 
gramming  of  the  system  matrix  and  is  programed  in  REAL  FUNCTION  S3. 
Since  the  data  XI,  X2  and  Y  are  already  provided  through  a  COMMON 
block  named  DATA3,  only  NPNTS1,  NPNTS2,  IP0WR1  (»  k-1)  and  IP0WR2 
(*  L-l)  have  to  be  passed  to  the  function  that  returns  the  value  of 


through  the  function  name  S3. 
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SUBROUTINE  MAT31  presets  a  submatrix  A^  in  the  following  way: 

i)  Preset  edge  section  elements  (using  REAL  FUNCTION  S3) 


4 

ii)  Copy  defined  elements  diagonally 

Q*;l« 
Ok 
•  •• 


°*j»  **jn  %a  a*) in 

%ya,yany-y>^ 

OtfJH 


<V-. 


rm 


Q  tt  I  O  Jr  pm 

T<+  «dgt  <^oc, 

hjU/  etiye  Stt*rO*l 


°*j» 

*•9*  aVL  <*kiU  — 


%  >M 

a*;l* 

?*/J!nr 
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%!! 

°Kn 

Qkjo 

Q*i  tl 

Q*;ll 

OifiU  • . . 

Gw, In 

3*; 31 

ak}3S  •  •  • 

•••  .  •  •«  _  *  •  »  »  •  « 
//✓A- 

ah,n 

a*j/h 

Q*;V 

a*at 

^k}U  •  •  * 

ty-A 

Qtyn 

air  hi 

QkjhS 

•  •  • 

Qkjnn 

This  subroutine  that  needs  NPNTS1,  NPNTS2,  M  and  IP0WR1  (•  k-1)  as 

* 

input  parameters,  returns  (SUBM1(7,7))  to  SUBROUTINE  MAT3,  which 
presets  the  entire  system  matrix  in  the  following  way 

i)  Preset  edge  section  submatrices  (using  SUBROUTINE  MAT31) 
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3.4.2.  Right  hand  aide  vector  B 

The  right  hand  aide  vector  B  will  be  divided  up  into  L 
subvectors  1.  class. 


Again,  the  asterisk  (*)  denotes  a  subvector  1.  class.  Now  the 

subvectors  1.  class  B*  can  be  written  as 

k 


2 

»  •  I 

2  EY„xir< 


,  L 


Finally,  the  elements  of  the  subvectors  1.  class  can  be  written  as 


2EY„xt'X2 


t-/ 

n 


To  compute  b^.^,  REAL  FUNCTION  S3  is  used.  Data  XI,  X2  and  Y  are 
available  through  a  COMMON  block  named  DATA3.  NPNTS1,  NPNTS2,  IP0WR1 
(■  k-1),  IP0WR2  (■  L-l)  and  IY  ("1)  have  to  be  passed  to  the  function 
and  the  value  of  b.  .  is  returned  through  the  function  name  S3. 
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>1*  ■ 


Using  FORTRAN  programming  language,  the  allocation  of  b 

k;L 

performed  by  two  stacked  DO  loops. 


can  be 


1-0 

DO  08  Ll-l.L.l 
IPOWR1-L1-1 
DO  08  Ml-l.M.l 
IPOWR2-M1-1 
1*1+1 

08  B(I)  -  S3 (NPNTS1 , NPNTS2 , IP0WR1 , IP0WR2 , 1) 
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3.5.  Software: 


The  software  to  compute  the  coefficients  for  a  three  dimensional 

approximation  is  described  in  Appendix  A  and  is  implemented  in  the 

TPL  HP-21MX  computer  system.  To  work  correctly  with  these  program 

modules,  the  user  has  to  conform  to  the  following  conventions: 

i)  Provide  the  data  in  three  arrays  (Type:  REAL)  of 
16  *  16  elements  through  a  COMMON  block  named 
DATA3 . 

COMMON  /  DATA3  /  XI,  X2,  Y 

REAL  XI  (16,16),  X2 (16,16) ,  Y(16,16) 

ii)  Dimension  an  array  (Type:  REAL)  of  7*7 
elements  to  contain  the  coefficients. 

REAL  COEF  (7,7) 

iii)  Define  the  parameters  NPNTS1,  NPNTS2,  L,  M  and 
IPRINT. 

NPNTS1  ...  if  of  XI  variations 
L  i  NPNTS1  <  16 

NPNTS2  ...  if  of  X2  variations 
M  <  NPNTS2  i  16 

L  ...  (desired  order  of  approximation 

polynominal  w.r.t.  XI)  +  1 

M  ...  (desired  order  of  approximation 

polynominal  w.r.t.  X2)  +  1 

IPRINT  . . .  controls  quantity  of  print  out 

2  ...  Print  system  matrix  and 

right  hand  side  vector 
before  and  after  Gauss 
Jordan  Elimination 

1  ...  Print  equation  system  after 

Gauss  Jordan  Elimination 
<0  ...  No  print  out 

>0  ...  Display  equation  (3.1.) 

with  the  actual  parameters 
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iv)  If  a  user  program  uses  subroutine  MAT3,  the  software 
modules  have  to  be  loaded  using  the  procedure  out¬ 
lined  in  section  2.4.  iv. 

If  all  these  requirements  are  met,  the  correct  call  for  the  subroutine 
is: 

CALL  MAT  (NPNT  S 1 , NPNT  S  2 , L , M , COEF , IPRINT ) 

Upon  completed  execution  of  this  approximation  routine  the  array 
COEF  contains  the  coefficients.  Externals  used  by  MAT3  are:  AB3, 
DATA3,  MAT31,  IEL3,  S3.  Under  no  circumstances  may  the  user  use 
any  of  these  names  for  modules  of  his  own  program. 


3.6.  Sample  User  Program 

3.6.1.  PTN^  Compiler  listing  of  sample  user  program 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 


>011 

0012 

0013 

111* 
0016 
0017 
08i« 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
•  026 

•  027 
0028 

m 

0032 

0033 

m 

•  037 

Soil 

•  040 
3041 
0042 
0043 
0044 
0045 
Cl  046 
0047 
0048 

•  049 

521° 

0051 

•  052 
0053 

Sols 

0056 
0057 
00S8 
0059 
0060 
061 


0064 

0065 

0066 

0067 


iiti 

6072 

0073 

0074 

0075 

0076 


FTN4 ,L 

PROGRAM  LAB  <3,99) 

*,Plot  calibration  point*  and  calculate  coefficients. 

COMMON  /  AFLD  /  PLOTR 
COMMON  /  DAT A3  /  Xi,X2,Y 
COMMON  /  AB3  /  A,B 

REAL  PLOTR (256) 

REAL  XI < 16 . 16) ,X2<16,16> ,Y<16,16) 

REAL  A<49,*9> ,6(49) 

REAL  YT  < 16 , 16 ) , R  < 16 ) , COEFF (7,7) 

INTEGER  IDCB(l44) ,IFILE(3> ,NOCR < 2) , ICLR < 3) ,YT1TLE<4> ,2T1TLE<4> 

DATA  PI  /3. 141593/ 

DATA  NOLF  /006S37B/ 

DATA  NOCR  /000033B , 040433B/ 

DATA  ICLR  /01SS24B, 01551 SB, 006537B/ 

C  FORMATS  SX^RT 

121  FORMAT  < •  Enter  raw  data  file  nane  "2A2) 

122  FORMAT  (3A2) 

123  FORMAT  < -  Select  1  ...  Pitch  angle 

Velocity  X  "2A2> 

124  FORMAT  < ■  Select  1  . „.  Pitch  angle 
8  -  <  ,  )  }"/ 

*  -  2  " ,  1H" , ' 

*nd  Delta  C  -  <  ,  )  >  *2A2> 

125  FORMAT  <"  Select  1  ...  Velocity  w/s  Ganna 
XX<  ,  >  3*/ 

*  »  2  ...  • , 1H" , 

*  C  X«X<  ,  )  J  "2A2> 


or  2  . 

u/e  Ganna  and  Beta 


u/s  Ganna  a 


and  Beta 


C  X« 


12640RW 


Enter  •  of 


160 

131 

132 

133 

134 

127 

128 

129 

130 

150 

151 

152 

153 

154 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
nax , 


Mach  Nunbers  and 
data 


<*  Enter  1  to  print 
</“  Calibrate  -  < 

</"  Calibrate  ■  < 

</"  Calibrate  X-X( 

</“  Calibrate  X-X< 

<24X*  *9X"  »9X"  *) 

<24X"  "9X"  *9X"  *> 

< 24X *  “9X"  "9X "X* ) 

<24X"  "9X“  "9X*X") 

values", 3<1X,F9 


"W) 
) "/) 
)"/) 
) "/) 


v/s  Ganna  and  Delta 
•  of  pitch  angle 
*2A2> 


155  FORMAT  * ( • 
X"  nax,  9 

156  FORMAT  <“ 


157 


*?6rma? 


<  "  mninun 

<"  naxinun  _ 

<"  of  user  coordin 
<"  Enter  1  to  redefine 
<“  Enter  nin,  nax, 

tartitions  ,  6 

nter  nin,  nax, 

partitions  ,  6 

enter  nin,  .nax. 


3)  > 


values", 3C1X.F9. 3) > 

bTes"$X"Xl"8X"X2"9X"Y“ 


artitions 
<*  Enter 


Xx.  •  partitions 


nin 


f6rma? 


158 

101  FORMAT 

102  FORMAT 
107  FORMAT 
105  FORMAT 

103*F5RMt 

X 

104  FORMAT 


4  x- 
i  X' 


) 

nax  ; 


these 

nin , 

’) 

nin , 

>) 

nin , 
n  in , 


data  *2A2) 


nax,  X  nin,  X  na 
nax,  X  nin,  X  na 


149 

611 

612 

613 

614 

615 

616 

617 

618 
602 


i  t  *  t  1 

( IX ,  3<F8 . 3 , IX, F8 . 3, 5X  > ,314) 

<"  Enter  order  M  of  XI  -  approxination  "2A2) 

<"  Enter  oredr  N  of  X2  -  approxination  “2A2> 

<"  Enter  IPR1NT  "2A2) 

<"  Select  1  ...  to  calculate  the  absolute  error V8X 
".4H"  ",9X"th#  rplative  error  "2A2) 

<*  Change  coefficients  ...  l"/ 

•  Change  axes  ...  2  "2A2> 

<  ".Enter  ALP MAX,  ALPHAY,  XO  and  YO"> 


FORMAT < " " <  <  3A2) ) > 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


<  "  1 " //// "  T It  va  mu*  d  ex  f  ex 
</16X,I2".  Mach  ntfnber"/) 

<6X "Ganna “6X "Be ta"7X"Ph i ■/) 

< 6X "Ganna “SX"Delta"7X "Phi"/) 

<  6X "Ganna ”6X "Be  ta "5X"X  vel*/) 

< 6X "Ganna “5X "De lta"5X“X  vel"/) 
< IX , 4  < 1 X , F9 . 3 ) ) 
</9X"Xl"8X"X2"9X"Y"> 

<  < 3X, 16  < 6X, 12) > ) 


f  r*  o  m  f  .1.  3.  «e 


rs  a 


wvb  rwni  irt  i  t  \  >  j.  w  t  wn  .  *  c.  /  /  f 

603  FORMAT  < IX , 12 , 16< IX , F7 , 3) /3< 3X , 16< IX ,F 7 . 3 ) / ) > 
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604  FORMAT  <///"  Absolut*  irror  at  «ach  point  C  error  »  (calculated 
Sualu*  -  measured  value)  )V) 

644  FORMAT  (///“  Relative  error  at  each  point  in  1  (  error  *  (calcul 

Sated  value  -  measured  val ue ) /measured  value  3"/> 

606  FORMAT  ( IX , 12 , 18( IX ,F6 . 1 )/3(3X , 18( IX ,F7 . 3)/) ) 

1101  FORMAT  <8H  Beta  > 

1102  FORMAT  (8H  Delta  ) 


1103  FORMAT  (8H  Phi  > 

1104  FORMAT  <8H  X  vel  ) 
LI  -  LOCLU(I) 

LO  ■  6 


read  rau  data  file. 


41  URITE  (LI,  121)  NOCR 
READ  (LI,  122)  1FILE 
URITE  (LlJ  149)  ICLR 

42  URITE  (LI,  123)  NOCR 
READ  (LI,  *)  IPX 

« “ ,o « 

IF  (IPX.EQ.l)  URITE  (LI,  124)  NOCR 
IF  (IPX.EQ.2)  URITE  (LI,  12S>  NOCR 
READ  (LI,  *)  IBD 
URITE  (Li!  149)  (ICLR, 1*1, 2,1) 

IF  (IBD.LT.l  .OR.  IBD. GT. 2)  CO  TO  43 

if  (IPX.E( !l  .AND?  IBD.EQ.i  )  URITE  (LI,  12 

IF  (IPX. £1.1  .AND.  IBD. EG. 2  )  URITE  (LI,  1 

IF  (IPX. EC. 2  .AND.  IBD. EG. 1  )  URITE  (LlJ  1 

IF  (IPX. EC  .2  .AND.  IBD.EQ.2  )  URITE  (LI,  l! 

CALL  OPEN  ( IDCB , IERR , XFILE , IOPTN , 0 ,28 ,144) 
IF  (IERR.LT.0)  60  T0’41 


*>  IPX 
149)  ICLR 

.1  .OR.  IPX #GT .2)  CO  TO  42 


IFLV  )W 

CALL  READF  < IDCB , IERR .X2.S12 .LEN ,5) 

if  (  ierr  ,lt.  o  )  Stop  0063 

RifB8F 

CALL  CLOSE  ( IDCB, IERR ,0) 

URITE  (LI,  126)  NOCR 

READ  (LI,  *)  NMACHjNPITCM 

URITE  (LI,  149)  ICLR 

IF  (NMACHfCT. 16)  GO  TO  44 
IF  (NPITCH.GT . 16)  GQ  TO  44 

88  IK  fcliWWtfi 

T(l,J)«YT(i,J> 


Cerrect  input  data  and  optional  output  of  data,  coefficient 


IF  (IBD.C^.2)  GO  TO  48 
DO  47  1*1 ,NNACH , 1 
DO  47  J-ifNPITCH,l 
X2< I , J)*X2( I , J)/X1(I , J) 

till*  <ti:  ,4Ji  "SSI  * 

WRITE  (LI,  149)  ICLR 
IF  ( IDUM.NE. 1 >  GO  TO  SS 

sra  M:n  mt\niA 

WRITE  (LO*  612)  I 

IF  (IPX.EQ.l  .AND.  1BD.E8.1)  WRITE  (LO,  613) 
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01S3 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

s  m  s 

0164  C 
0165  C 
0166  C 
0167  C 
0168 
0169 
0170 
0171 
0172 
0173 
0174 
0175 
0176 
0177 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
0185 


IF  (IPX.EQ.l  .AND.  IBD . EQ . 2 )  WRITE  (LO,  614) 

IF  < IPX . EQ . 2  .AND.  IBD.EQ.l)  WRITE  < LO  615) 

IF  < IPX . EQ . 2  .AND.  IBD.EQ.2)  WRITE  (LO,  616) 

DO  S4  J“1 ; NP ITCH , 1 

54  WRITE  (LO,  617)  XI ( I , J ) ,X2( I , J > , Y ( I , J ) 

WRITE  (LO,  618) 


.Calculate  nininun  and  naxinun  values  for  all  coordinates  and 
optional  redefinition  of  itiniitun  and  naxinun  values  by  user. 


55  XUHAX-Xl(l,i) 
XUMIN-X 1(1,1) 
YUHAX-X2(i,i> 
YUHIN«X2( 1 . 1 ) 
ZUHAX*Y (1,1) 
ZUHIN«Y(1,1> 


DO  45  I-1,NHACH,1 

do  45  j«i;npitch,i 

IF  (X1(I,J) .GT.XUH# 


IF  (X1(I,J) .CT.XUMAX)  XUMAX-Xi(I,J) 

IF  (X1(I  J) .LT.XUMIN)  XUHIN-X1(I,J) 

IF  (X2(I,J) .GT.YUHAX)  YUHAX-X2( I , J) 

IF  (X2<i;J) .LT.YUHIN)  YUMIN*X2( I, J) 

IF  (Y(I,J) .GT.ZUHAX)  ZUHAX«Y(I,J) 

IF  (Y(IjJ).LT. ZUHIN )  ZUHIN-Y(I  J) 

45  CONTINUE 

IF  (IPX.EQ.l  .AND.  IBD.EQ.l)  WRITE  (LI,  127) 
IF  (IPX.EQ.l  .AND.  IBD.EQ.2)  WRITE  (LI  128) 


IF  (IPX.EQ.l  .AND.  IBD.EQ.2)  WRITE  (Li;  128) 

IF  < IPX . £0.2  .AND.  IBD.EQ.l)  WRITE  (LI,  129) 

<Lt;  iSo> 

WRITE  (LI,  151)  XUHAXjYUHAX, ZUHAX 
WRITE  (LI,  152) 

INCX  »  4 
INCY  -  4 
INCZ  «  4 

WRITE  (LI,  153)  NOCR 
READ  (LI,  *)  I 
WRITE  (LI.  149)  1CLR 
IF  (I.NE.l)  GO  TO  49 

IF  (IPX.EQ.l  .AND.  IBD.EQ.l)  WRITE  (LI,  1§4> 

IF  (IPX.EQ.l  .AND.  IBD.EQ.2)  WRITE  (Llj  155) 

IF  ( IPX . EQ . 2  .AND.  IBD.EQ.l)  WRITE  (LI,  156) 

IF  ( IPX . EQ . 2  .AND.  IBD.EQ.2)  WRITE  (LI,  157) 

READ  (LI,  *)  XUHIN,XUHAX,YUHIN,YUHAX, ZUHIN, ZUHAX, INCX, INCY, INCZ 
WRITE  (LI,  158)  XUHIN,XUHAX,YUHIN,YUHAX, ZUHIN, ZUHAX'INCX, INCY, INCZ 


.AND.  IBD.EQ.l)  WRITE 
.AND.  IBD.EQ.2)  WRITE 
6)  XUHIN,YUHiN, ZUHIN 

1)  xuhax;yuhax, ZUHAX 

2) 


2)  WRITE  (Li;  128) 


020 3 
0204 
0205 

m  g 

0208  C 
0209  C 


Initialize  plotter:  define  plotter  area)  calculate  scalinq 
coefficients  XA,  X8,  YA,  YB,  ZA  and  ZB. 


49  CALL  INITG  (13) 
XPHIN  -  0 . 
XPHAX  -  10. 
YPHIN  ■  0. 
YPHAX  »  10. 
ZPMIN  -  0. 


ZPHAX 

ALPHAX 

Mi 

xx8 

XL 

YA 

YB 

YL 

ZA 


<<^S!S«xuHAx:^^«xuHiNl//^ra:SB^^ 

(XPHAX  -XPHIN  ) 

(YPHAX  -YPHIN  )/( YUHAX— YUHIN) 

( YPHIN«YUMAX-YPHAX»YUHIN  > / ( YUHAX-YUHIN ) 
(YPHAX  -YPHIN  ) 

(ZPHAX  -ZPHIN  )/ (ZUHAX-ZUHIN) 
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ZB  »  < ZPMIN*ZUMAX-ZPHAXSZUMIN> /( ZUMAX-ZUMIN) 

ZL  ■  (ZPMAX  -ZPMIN  ) 

XO  ■  8. 

YO  »  12. 

BS  PL0TR(61)  -  < ALP HAX*P 1 ) / 1 8 Q . 

PLOTS  (62)  -  ( ALPHAYSPl ) /ISO . 

PLOTS  <63 )  *  ( ALPHAZtP I )/180 . 

PLOTR  <  64)  -  XO 

PLOTS  <  6S )  -  YO 

PLOTR  <66>  -  XA 

PLOTS < 67)  -  XB 

PLOTS (68)  -  YA 

PLOTS (69)  -  YB 

PLOTS < 70)  -  ZA 

PLOTR (71)  »  ZB 

CO  TO  (86,87)  IBB 

86  CALL  CODE 

WRITE  ( YTITLE ,1101) 

GO  TO  88 

87  CALL  CODE 

WRITE  (YTITLE ,1102) 

88  CO  TO  (301,302)  IPX 

301  CALL  CODE 

WRITE  (ZTITLE, 1103) 

GO  TO  303 

302  CALL  CODE 

WRITE  (ZTITLE, 1104) 

303  CALL  SETSM  (113.1.) 

CALL  AXIS  ( XO , YO , XL , ALPHAX , 6H  Ganna,  e.XUMIN.XUMAX , 4HF6 . 2 , 6, INCX ) 
CALL  AXIS  ( XO , YO , YL, ALPHA Y^  YTITLE , -8, y6mIN,YOMAX , 4^F6 . 2,6 > lACY ) 
CALL  AXIS  ( XO , YO , ZL , ALP HAZ J ZT I TLE ,  -8 , ZUM IN, ZUMAX \ 4HF 6.2,6,INC2> 


(ALPHAX*PI)/180 . 

<  ALPHAYSPI ) /180 . 

( ALPHAZSP I >/180 . 

XO 

YO 

XA 

XB 

YA 

YB 

ZA 

ZB 


Plot  neasured  calibration  surface . 


CALL  SETSM  (113,2.) 

DO  n  1«1,NMACH,1 
DO  01  J-i,NPITCH,l 

CALL  THRTW  (XPLOT , YPLOT .XI ( I , J > ,X2( I , J ) ,Y( I , J > ) 
IF  (  J  .EQ.  1  )  CALL  PLOT  tXPCOT, YPLOT, 2) 

oi  hWjiU  ^-  1  >  CALL  Pl0T  ‘’"“-ot, yplot, 3> 

DO  02  J»1,NPITCH,1 
DO  02  I-l,NMACH.i 

W-V  ■«teTi>f!>)1  itttH'.MWtt  : 1 

IF  (  I  .CT.  1  )  CALL  PLOT  (XPLOT, YPLOT, 3) 

02  CONTINUE 


Calculate  calibration  surface  coefficients. 


91  CALL  SEtSM  (113,0.) 

WRITE  (LI,  101)  NOCR 
READ  (LI,  *>  MORDER 
WRITE  (LI.  149)  ICLR 
IF  (MORDER . GT.6)  GO  TO  91 


?a  till1  (ti; 


10?1  » 
tm  ttii  m  m 

BKtfc  14?1  wr 


IF  (NQRDER .CT.6)  CO  TO  92 

M*MORDER*l 

N*NORDER+t 

CALL  MAT 3  (NMACH,NPITCH,N,N,COEFF,IPRINT) 
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Overwrite  data  array  with  calculated  data,  based  an  the 
just  obtained  coefficients. 


WRITE  (LI,  105)  NOCR 

READ  (LI,  *)  IAR 

WRITE  (Li!  149)  (ICLR, 1-1,3) 

IF  (IAR.LT.l  .OR.  IAR.GT.2)  GO  TO  411 
IF  (IAR.EQ.l)  WRITE  (LO,  604) 

IF  (IAR.EQ.2)  WRITE  (LO.  644) 

WRITE  (LO,  602)  ( J , J-i , NPITCH , 1 ) 

DO  20  I-1,NNACH,1 
DO  1?  J-l, NP ITCH, 1 
SUM-0 . 

DO  96  I 1-1 , M , 1 
IEXPl-Il-i 

IF  (IEXP1.EQ.0)  GO  TO  401 
X1EXP-XK  i  J)**IEXP1 
GO  TO  402 
X1EXP-1 . 

CONTINUE 

DO  96  Ji-1 ,N, 1 

IEXP2-J1-1 

IF  (IEXP2.EQ.0)  GO  TO  403 
X2EXP-X2  ( 1 , 1 )  MIEXP2 
CD  TO  404 
X2EXP-1 . 

CONTINUE 

SUM-SUM+COEFF ( I 1 . J 1 ) SX1EXP SX2EXP 


SUH-SUM+COEFF ( I 1 . J 1 ) 

HjJSMlMW  iA* 


GO  TO  19 

R  < J )«< (SUM— Y (X,J))/Y<I,J))*100. 

Y ( I , JJ-SUM 

WRltE  (6,603)  1  ,  ( R  ( J  > ,J—1, NPITCH , 1 ) 
WRITE  (6,602)  (J, J-l , NPITCH, 1 ) 


Plot  calculated  calibration  surface. 


CALL  SETSM  (113,3.) 

DO  81  1-1 ,NMACH . 1 
DO  81  J-i,NPITCH, 1 

SPViSTij  ,Stt*»Z5V,T*KW:  >  -T<  1 

8,  {WMl-11  C*LL  ,L0'  <"M-0T>L0T,3, 

DO  82  J-l, NP ITCH ,1 
DO  82  I-l'NMACH,l 

5P*-YiT8ST?  >  *  £ShiT»t5V°  <  ;  iicSf  <  T  *  * . ■»*  > 

IF  (I.GT.i)  CALL  PLOT  (XPL0T,YPL0T,3) 

82  CONTINUE 


Redefine  Y . 


i>o  6o£  i-i ,  nma£h ,  1 
DO  802  J-l 'NPITCH 
802  Y(I,J)-YT(I,J) 
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0381 

0382 

0383 

0384 

0385 

0386 

038V 

0388 

°oll9o 

0391 

0392 

0393 

0394 

039S 

0396 


C 

C 

C 

C 


N«Xt  ST«P  ? 


WRITE  _ 

READ  (LI,  *>  _ 

WRITE  (Li!  149)  < ICLR ,1"1 ,2,1) 
IF  (IDUM.EQ.l)  GO  TO  91 


(LI, 103)  NOCR 
(LI,  *)  IBUM 


IF  (IDUM.EQ.2)  GO  TO  89 
CALL  STOPG 
STOP  077 

89  WRITE  (LI, 104) 

READ  (LI,  *)■  ALP H AX, ALPHA Y,XO,YO 
WRITE  (Li'  149)  ICLR 
GO  TO  85 
END 


FTN4  COMPILER:  HP92060-16092  REV.  1926  (790430) 

t*  NO  WARNINGS  t*  NO  ERRORS  f*  PROGRAM  •  0446S  COMMON  «  00000 
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•  Load  mao 

points  and  calculate  coefficients 


to  load  these  program  modules,  enter  (from  LOADR):  MS , itTPLBL 
MATS3  mil  ?n"f25rox-','Qtion  7  DATA3  7 

MAT31  ^74?I  2721?  IbIaSSE 2 ! t"a! *on  arr*n<}«  »y«te*  matrix  and  vector . 
S3  27652  30070  IdIaSS^J arranqe  submatrix  i.  class. 
llz  5otll  32i5S  MSPSHiSaM  /  AB3  ,  sunmations. 


.  ur aw  and  label  a 

ol  array  Ai  initialize  plot 
nove  pen  to  a  defined  po 
chanqe  plotter  mo 
tcrMinate  qraf 
coordinates  to  plotter  un 

vary  size  of  the  p 


*  mm 

3  790215 
3  781229 

3  780S26 
3  790202 
3  790223 
1 
1 


3.6.3.  Reaulta 

i)  printed  output 
Thaaa  d  a  t  n  f'r  om 


Ganna 


- .  810 
-.519 

-:?1S 

-.103 

-.059 


Ganna 


-.781 
-  .506 
-.280 
-.178 
-.092 
-.048 
-.004 
.  032 
.061 
.130 


Ganna 


-.787 

-.510 

-.276 

-.172 

-.085 

-.0S0 

-:°o$* 

.066 

.131 

•  m 


Ganna 


::??! 

-.276 

-.043 

-.003 

.93? 

.064 


1.  Mach  nonber 
Beta  Phi 


35.000 
SS. 000 


-5.000 

-2.500 


i:iu 


5.000 

10.000 

£:in 

35.000 


2.  Mach  nunber 


-3S.000 
-25.000 
-IS. 000 

-lg.aoo 
-5. 000 

“o’.ooo 

i.so°o°o 
10.000 
IS. 000 

li.000 


3.  Mach  nunber 


-I: loo 

1:181 

5.000 

10.000 

hi.m 

35.000 


n 


file  e<«»L_  X 


-35.000 

-25.000 


-15.000 

-10.000 


4.  Mach  nunber 


--iim 

-15.000 

"-§:!!! 

-2.500 

0.000 


!:??! 


10.000 

15.000 

25.000 

35.000 


Ganna 

Bata 

Phi 

-.795 

-.045 

-35.000 

-.521 

-.  045 

-2S. 000 

-.282 

-.045 

-15.000 

-.175 

-.048 

-10.000 

-.082 

-.050 

-5.000 

-.040 

-.050 

-2.500 

-.010 

.034 

=:»» 

.067 

-.049 

5.000 

.134 

.198 

.350 

-.050 

-.053 

10.000 

15.000 

-.  055 

25.000 

.593 

-.050 

35.000 

6.  Mach 

nunbar 

Ganna 

Bata 

Phi 

-.803 

- .  06S 
- .  06S 

-35.000 

-.528 

-25.000 

-.£64 

-.065 

-15.000 

-.173 

-.068 

-10.000 

-.085 

-.  071 

-S. 000 

-.032 

-.072 

-2.500 

-.004 

-.071 

0.000 

.  024 

-.071 

2.500 

.058 

-.071 

S.  000 

.131 

-.071 

10.000 

.199 

-.075 

IS. 000 

.349 

-.078 

25.000 
3S. 000 

.603 

-.073 

XI 

X2 

Y 

45 


272E-M  .777E-IS  ,»33E.*1  -.312t.lt  .MSt-l 


*  «*  n  u»  s  rs  m 


«  <  i>*  ki 


•  4  CD  *•»««•£ 

f*  *  N  «•  ft  «•  4  N  «• 


—  ri»  i/»  —  r*i  »  u*  ■-  r*i 


Uj  Ul  Hi  IK  UJ  K  111  IK  uj  U  J 


Uik^wJUJuJWuKUUJUiuJ  M  jJ  IK  IK  iK  4l  IK  IK 

n3b<«im»i<i3«nn  K«4ui>»ii 


■«  p*.  »  • 


pv  ©  >  • 


IKKKIKKWIKKKIKIKIKW  K  K  Ul  K  IK  K 


<  N  *  M  *  Ifl  t* 


tii  s 


I  ♦  I  ♦  I  ♦  I  I  I  I  !♦!♦«♦(  I 


ssawfisisrirssijRsssicgsigssas: 


4. 


FOUR  DIMENSIONAL  APPROXIMATION 


4.1.  Problem: 


*2 


A  data  set  of  NPNTS1  *  NPNTS2  *  NPNTS3  data  points  is  given, 
where  Y  depends  on  the  parameters  XI,  X2  and  X3.  The  data  pattern  is 
to  be  approximated  by  a  function  Y  -  f(Xl,X2,X3),  so  that  the  error 
between  data  points  and  analytically  determined  points  is  lower. 
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4.2.  Approach 


Again,  as  in  the  two  and  three  dimensional  case,  a  polynomial 
method  is  used.  For 


Y-  f  CXI ,  xz,xi) 

we  assume 

Y-  fc,  *  cnx3*  - 

N~f 

*  <*,»'  - 

.  6a.  tf"  7-  >3  * 

•  0  ♦ 

A/-*  M-l  ■% 

/V-/ 

7* 

...‘C^J.Xl  * 

•  ♦  • 

-*  ifdz  -»  C,MJ  x/-# 

/v-/_  .  7  . ,  .  , 

7-^  7Ar/ 

-»  ccui  f  Cuuv3  i  ...  *  Cu.nX3  7*  *-3  ■* 


*  [  Ct -  + 


'Ln  L*i 


or 


#-/,  ***?„,  L‘ 
*...'CLMMXS  J-XZ  fX! 

l  n  H  K-i  j-i  , 

Y.  E  { E  f  E  C-J.X3  }*2J  }  xi 


(4.1.) 
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This  data  arrangement  which  discriminates  between  data  points  taken  at 
constant  XI,  X2  and  X3,  respectively  is  needed  in  the  application  of 
this  method  to  the  description  of  pneumatic-velocity  probe  characteristics. 
The  various  surfaces  in  the  above  sketch  for  example  can  be  seen  as  data 
taken  at  different  Mach  number  (XI).  At  each  constant  Mach  number 
(XI  “  constant)  yaw  angle  (X2)  and  pitch  angle  (X3)  are  varied,  and  the 
probe  is  measured  to  be  Y.  Even  though  probes  are  usually  balanced  to 
the  average  flow  yaw  angle  in  some  cases  this  is  not  possible  and  the 
four  dimensional  approximation  is  then  needed.  A  similar  method  for 
describing  probe  calibration  data  is  used  at  the  institute  for  Jet 
Engines  and  Turbomachines  at  the  Technical  University  Aachen  (Ref.  2). 

The  least  squares  criterion  is  then 

mmci  nwhu,  mjwim 

t  fCxtrH  /  x2rrl  /  rsl)  ~  Kst  ^  2  ) 

where  the  indices  rst  denote  data  points.  Using  equation  (4.1.),  R 
becomes 

jpiftti  Hfm fW  M*l 

rj  ^  C'»4  *  Ciin  *3r*t+ 

*  (c*  ■*  c-«  x3y  •••♦  cu*  * 

•  •  % 

*  4  —  4  ^  ihn  ^  *  X&nL  * 

cin  4  C«  +C2»H  X^nt  * 

'(Cm  ♦  *.  .  4Czxn  *2N~'  + 

•  •• 

4(‘w.  +  Ct»**J  tim’d"  )'x£t'rXlr* 
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+  rt*2«t4  •••  4  C^yl 


Ccy,  ♦  '*•  4cl 


*3  )*^2-*  + 


The  term  In  the  [  ]  -  bracket  shall  be  called  B 


IMSt  AWWflf 


*-E{Lf2  «TH 

r*/  <»< 


‘  'O#  >*  ',v; 


4.3.  Solution: 


To  minimize  the  error, 
to  the  coefficients  CJA,  and 
thus: 


R  is  partially  differentiated  wit!  respect 
the  partial  derivatives  are  set  to  zero. 


™  l0 


/*  /  ■■■  i  ;  jc  /  ..  .  M  j  k*  /  ...  /V 


34* 


nmrsi  N*m r/ 


Cf£{  E  2eMu-o 

rv  •»■/  *»/ 


Assuming,  that  the  summations  extend  over  all  data  points,  ZEE  should 

«*!»*'  wmi 

be  understood  to  be  Xj  Jj  VJ  .  Now  even  a  very  large  sheet 

r*/  *./  i., 

of  paper  can't  take  the  system  matrix,  unless  we  skip  to  submatrix 
notation. 


fa  II  ^ 12 

fall 

A* 

fa it 

<c; 

*1 

s* 

s* 

•  •  • 

K  ... 

K 

m  •  + 

< 

•  •• 

< 

(4.3.) 


The  asterisk  indicates  a  submatrix  1.  class.  On  the  following  pages 
the  submatrices  and  subvectors  1.  class  are  listed. 
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/*>nm  *  burnt*  m»v?j 

EEEw* 

EEE*C 

E££w„, 

I E  Ew„* 

•  •• 

e  r 

ECE< 

eee«; 

E  E  E*^,**?,* 

EE  EWb,**^, 

•i« 

E  EEw^.xC 

EEE*C 

EEE*£ 

eeeC 

EEEw^-xC 

EEEtt,*-xl,* 

«•« 

ECExi^xC* 

EEEw^ 

CEE«*<' 

EE£x£, 

EEExi*-m„, 

EE  £«?„,?  *£' 

ELE*.*^ 

EEE 

E  ££«„.*£ 

EEEw*-*^ 

EEEm^xjJI, 

i  •" 

EEExiV’^w 

m*v< 

ErE*z_**j^ 

E£E*v*C 

EEE<- 

. 

eee*i;. 

.  ••• 

EEE<-<* 

i 

1 

•  •• 

- 

i 

EEE< 

£  £  E  *C*xi* 

CEExC-C 

ECExC 

EEExj,*-  xjm 

EEEtt^xC 

EEE<-*4, 

£££»*•»& 

EEE*C-*C 

| 

EEE<-«U 

E  E 

££&£.«£ 

EEE«' 

E  E  £  --• 

EEExOC'1 

EEE<-*C 

EEE«**«« 

EEE*z£i*x£l 

i 

i 


I 


*£•««» 


£EE*C-*C 


LlZ*C-*Cl 


FlZxlC 
EEC  xC*i* 

EE  ErtlnT*^* 


EEExC-«m 

EECtt^.xj^ 


E  E 


EE&a^*-xC 

££DC:«; 

EEE<?lc 


EEEXi„ 

EEExi*-** 

CUE  xLr*J»*  ••• 

CEE  *U-K 

EEL 

EEE*w*& 

EE  Exi^xi^, 

EEC *i*t 

EE  E  xi^«,^-«nt 
EEExinix2^-xi^l  ... 

CCEx^  xi,.  *£  j 
EECxijx^.w,*  ; 

CEE  xi,**#^ 

EES^-xj^ 

E  E  E  x  xjul 

EE  Ex'^xi^x?^ 

EEC*1*  xiw->^  , 

EEExu-x:* 

EEExu-^-xi^ 

E  E  Exu-xWh* 
ErE>..-*i^-xC 

EEExi— ‘XZ^-xJ,,, 
LrExi-xi.r*C 

t 

EEExi„-«i, 

EEE*'H.<.-xi,ri 

EEEx^xC^  ... 
EEExi„,x£,x/^  ... 

EEEvrc*-'# 

EECx^-xC-xC 

EEExu-xi^xC' 

EEExi,*  xi^-xj’ 

E  E  Ex',*-  X^nt'Xj^ 

ESE  *'„■« 

EEExin(-x2H(-x.?n'  ... 

eecxv»4-<‘ 

•  •  • 

*  —  "  ’ 

r 

_ .  i 

EEEx^.xC 

EEExvC-xj* 

EEEv^'-w^ 

E££xim.»£.«^ 

CEE  xi  ‘XiM 
EEEv>C< 

EE  Exiw-xi,* 
EEExi^-xii  xsM 

EE  rfutCxt-i  ... 

£  EEx^jfi^x/^ 

£EExiM,-x£xfcI" 

EE  Exi^.xz^.  x^ 

EEExu-<;-w^ 

EEExvxC-xC 

EEExw*vC 

EE  Exi„  *£„,  xj,^ 

EEE*'*  x£  xj;  ... 

E  E  Exi^-xi^-x/rt, 

I 


I 


I 


C  E 

... 

ccr  x^-xi,,- 
EECxiHt-xi,,-x/,'t 

... 

EEE»*-< 

E  EE  X3— 

EEL  x^-xj^-  xj„ 

EEEvim-*C'-x4i 

E  E  Ex^-xi^-tf,* 

£  £*'■<  x3„i 

EECx^xz^x-C1 

EZLfi^nZ-tC 

EEEWjxi^-k^ 

EEn'^-xl^ -xim 

CE x^xi^  xi,^  ... 

CExux4,«t, 

EEEVtt',.*# 

EECxi*.*^*^ 

i 

i 

( 

EEExw-xz- 
E  E  Eyi^-xz^-xi^ 

EEExi^<-xj,t  .. 

E  E  Exit<  •  xz^’XJj^ 

EEExu-<  x£ 

E  E  ^n* 

EE>i«r«ii-w£  ... 

EEE>W*&»C 

E£Exi„-xzl  -*C 

EEE*'*«  - 

EEE»U«;*C 

h  !  ' 

j 

h  '  i 

i 

:rxu  xi;,  «.i  ... 

... 

E  I  £»«<•*£■ 
EEExirt-«”.rC 

1 

! 

"EEE^.xC 

CEE^C’-*}* 

EECx^'*^*^ 

E  £  Exi^xz^-wJ^ 

nrxW.xC«r 

lEEx^-xC*  < 

:cv«i<  ... 

EEEx,„.xi:.xC* 

EEExu-*C< 

E  CEXl^  xz^, 

EEDU-xC^C 

I 


S?' 


a 


V.- 


EEE£ 

EE  ZxC-»j„ 

EEE«"-*'« 

EE&C-«fc 

EEExC'-xC 

ZZExC-xi* 

EE  E  Jrt^-AZ^i 
EEExi^x^, 

EEExC-xZ^xj,,,  ... 

E  E  Ex^t-xz,*  -w^, 

E  E  Exi,rtxzw  x/m 

EEExCX*C 

EEDw«C 

EE&CX 

EEE»C«r 

EEEv«.,*xC 

EEExCv*^  - 

EEEOV'C* 

EEE*£x 

EEExCxX 

E  EE  x£ 

EEExCw^^i 

E  ££<•«*■*- 
EEExtXX 

EE  £<•«*, 
EEE/C-«-*4* 

EEExi'Xx  ... 

E  E  Exi„,  xzJrtXJ,^ 

EEE«-«r 

EEE<C.  <i  <t 

j 

EEECVC 

EEExiS^-xi; 

1 

EEExCXxi,7 

eeexC-*4*C 

eee<|-»4-*C.  ..." 

EE  £«*<*#* 

••* 

-■■-■»-  - 

.  -  • 

—  -  -  •  . 

| 

EEE<-*C 
ZEExC  xC*1,* 

E  C Exi*  xi**  *JTgt 

EES<K-*A  ... 

EEDCxC  -xC 

EEE*CC-<C 

EEEC*C 

E  E  £«•«,,, 

EEExi^  xz^  *J„, 
EEE<<x4  ... 

EESC-*C-»C  1 

EEE<«-x3- 

EEE^xCC 

EEExi^-u^xj^  ... 

EEE>C-*C-*tf* 

EEExC-xC-xC 

E  E  Ex'-x2Vj1  •  xJ— 

EEZ<<<‘ 

<<'•«*■  x*,,,  ... 

S  E  TjXiff<  *2^ 

... 

EEE*C,;-*C< 

•  ## 

... 

ECE^A*-*#* 

... 

eee«  «r 

_.  . 

<;x&x&  ... 

SEECxii,  <4 

*•• 

*'r**^w* 

EEExC^xiiT* 

■S-»C-*JU  - 

T.ZZ«m  *C<’ 

... 

L'EE>C«>3i 

C*C«I  - 

EEr<;  *2;  <-» 

222*W*C 

E22*C*0»** 

2EZx£  *C«W,  ... 
EE£x£»j£*4  - 

ZZEx&xC-xC 

EEExC-«C*Ji 

EEExCC*C 

£££«•<  ".* 

EEExi;-*-.^ 

E  EExim.x2£ 
EEE>C  «^  x3^ 

EEExi^xi^xj,^  — 

EES^-x^-xi^  ... 

EESO-xC1 

E  E  £  ^1*  *xz^  •  x3„i 

EEExi^-ttJJ,  xi^ 

E£Exir(  *V  xj,*  ••• 

TlLDC-iC-nS4 

- . - 

- - 

EEEx>:<‘ 

'  EEExfr«5X. 

EEExux^.xj^ 

EEE<  «£*•*&  .- 

MEC*C-«C 

EEE^-xC-w; 

EEExC-iC-aS 

EEExC«i^xC 

EEExi^'-xC1 

EEExC  : 

EEExii,  x^, 

•  •  • 

EE.Exii.xj” 

EEEX-xz* 

£EExiitxi,ilxj„,  ... 

EEExi, *2,t 

CEExiJeXj* 

CEExi*X 

#«# 

EEEXX 

EEExii  xi^  xi^ 

EEExii,^)®^  ... 

EEEX  x.„ 

i 

EEExC<' 

EEE>i*x£ 

EEEXX*"‘ 

EEE*i,iiXilltxj” 

EEEX>z„,X  ... 

'  EEExti,^ 

CEEx'ix^i 

E  E  Ex.^.  x2^  ■  x3^ 

««• 

EEijXi^ 

EEExi^x^ 

EEExCX.x4,  - 

EE  EX  xi\ 

EEEx^xz^  x?* 

EEZX  w*  «* 

... 

EEExi  i,>anlX 

EECX*^*^ 

EEExii.  tti,  X 

EEEX*2! 

EEE*C«n«XC 

EEE*|^,x.zHlxJlll 

... 

SEExii.xi^x^ 

EE2X«T 

EEEX<,X  ••• 

EEExiX 

;  ; 

9  99 

1 

1 

- - 

'  ! 

:  i 

i . : 

. 

_ 

i 

I 

.  .  .. 

EEE*C*C 

E  E  EX  x£  xi* 

... 

EEExC<-xC 

SEExCC 

EEExCXw*  - 

EE  EX  xi, 

EEZ>CxC«- 

EEExii,xCx^ 

•  M 

EEExCxOX 

EEExii.-tti.xJ^, 

EEExiiXX 

EEExiJgxii 

EEExCXxC 

EEExC-xS-xC 

••a 

EEEXXXT 

EEEX»£tf* 

EEExi„xitil 

EEExi^xi, 

tl 


( 


... 

E£Ex!,rt  x2m 

»'*«.*■  wit  ... 

EEE*C«*«,S 

•  •  • 

*c  ... 

CE2X  x2‘«*3*'2 

EEExiJX«C 

. :  . . 

"*<»  *4  ... 

EEE«i-n^«,J 

•  •  • 

eee^-:  <«r 

1  - - 

eee*c< 

E  E  Cx^  x£  tivt 

EEE*C.*C*C 

ZLLx'+Kt «* 

EEC Jcfe-KSaU 

...  EEExiixijxi^ 

EEDCxC»C 

EEE<t<'< 

...  EEExwx2*  x./^ 

EEE*C-*C _ 

EEE/i;txC«rt 

- EEExii.xC'C* 

DSEx^xC** 

••• 

EEExi*  *Z*  xj^ 

...  EEEjm^jo,,,  is^t 

EEI>C,-Mi-xC 

1 

EBEWirWW 

i 

!  EEExi^xC  xj^*4 

«»*  ... 

EEExi^  xiw  xjM 

EEE*i,-xC 

EEExC-xC-xirti  - 

EEE^xC^r 

... 

EEExi^  *C«it 

•  •  • 

EEExi*  x<5*  x&  ... 

EEE<,*C*C 

... 

EEE^-aiwi** 

SEEWUOC 

E  C  Exi^xz^  x^  ... 

EE£x£x£4-x£J 

£  2L_ 


EE2xi“ 

£££*£.  xsM 

EEExiJ^xi* 

EEExi^xj^ 

EEaC»r 
EEExtf  < 

EEExij^xi* 

ECEmS1^^ 

zjEEx^’xi^i  xj* 
EEEx^xi^xzit  ... 

EEEx^xi. 

CEExCX 

EE&C»C 

EE&C-xC 

EEExtfxC 

•  •  • 

EEExC  xi« xZ^, 

EEExCTx^-x^  ... 

EEEC^ 

EIExi^'xi^ 

EEEJC-’Vtf* 

EEEx^1  xz^  x^ 
EEExiy-w^-ifli  ... 

E22*£  xz^xC 
EEEC**-«C 

zmc,'< 

EEExi^xz^  xi^ 

ZEExi^'xz^xJ^  ... 
££&#<,*&  ... 

EEExC2*^ 

E  E  £xi*  «?* 

EEE<X< 

EEEx^xz^  xj^  ... 

EEEx^x^xC* 

EE  &<£<■«£ 

EEExi“xi^x?^  ... 

2EE)C'<. 1 

•  •• 

-----  ■  - 

-  - 

- 

EEEx£*C 

EEExffxC** 

1 

EEExffxC*-  ... 

ee&c?xc-*4  ... 

EE&C«S'«C 

EEEx£*xz£  «£ 

EEExtfw; 

EEE>C*0*U 

E  EDxi^  x£  x4, 
EEExi^xz^xii  ... 

EEExC?-xC: 

EEExlJf-  xiw  j 

EEDff<*C 

SEE***  x2^  x^,,  ... 

E£E*ff<;  «£* 

EEExi“<,x£ 

EEEx£J<lxC  ... 

EEExi“  x£  ; 

LZlY^-xi^ 


The  double  asterisk  Indicates  a  submatrix  2.  class.  All  submatrices 
2.  class  on  the  same  diagonal  bands  2.  order  of  the  aubmatrlces  1.  class 


(k  -  1,...,2L)  are  Identical  and  therefore  will  be  renamed  from 
**  ** 

“Scjitj  t0  “St;! 


V  - 

A** 

'V,' 

• 

A4* 

*  2,1 

*  A** 

r'k;!.Z 

.  as 

*lw  ■ 

A*4 

■  A** 

/'k,2.Z 

A  44 

*  Alt;  u 

»  •  • 

A**  * 

A  ** 

Ak,M./ 

.  A " 

rk,tH,2 

*** 

*  >444 

•  A*-’/. 

K*« 

A** 

’A*;#.* 

M** 

/h  ^ 

'  <x» 

K 

*  '  *  <r;  +.  n-t 

A*4 

•  me 

z*‘ 

44 

*  ^  /,* 

where  A,  can 

I  Z2E« 


can  be  written  as 


ZE&C<h, 


«4 

W  ..,2 


m  ac«:*c 

ll  BOaK-*C 


km  ...  ,  2L-/ ;(•  !,  ..v  2M-! 


SELX.XC 


r  / 


Note,  that  this  submatrix  2.  class  Is  not  only  symmetrical,  but  also  the 
elements  on  each  diagonal  line  2.  order  are  Identical.  The  elements 

j  o£  the  submatrix  2.  class  A^  therefore  are  renamed  to  a,  ; 

9  •**  k  ;  L  f  IS 

L  Specifies  the  place  of  the  submatrix  2.  class  in  the  submatrix 
1.  class  A^,  k  specifies  the  place  of  the  submatrix  1.  class  a£  in  the 
system  matrix  A  and  i,j  specifies  the  place  of  the  element  a.  in 


the  submatrix  2.  class 


** 

V,i 


*  R.  . 

4/ 

ct,  , 

-ct  -  a. 

hlil 

k;l;2,l  *il>  U 

• 

a*,if2 

««  • 

ak;l,2,t  '  ak,l,t.2S 

*  %/;/</  V;W/  *  ** 

akyl;M*f 

'%i,W  ■  ■  • 

'  *  a*;02.A( 

s  a*,4  • 

k^ttH 

-her. 

m  can  be  vrltten  as 

a.  .  -“2 {'SiE'xiX-xg.xsZlt 

k;L)m  »*«i  fc»i  W  Wt 

>J  -  ,-2£-'  ;  'y  ....  y3M-l  J  hn.  ty 


Kjlf  H,N 


(4.4.) 


.../My 
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This  operation  is  programmed  in  REAL  FUNCTION  S4  (NPNTS1 .NPNTS2 , NPNTS3 , 
IPOWR1 , IPOWR2 , IP0WR3 , IT) .  The  data  XI,  X2,  X3  and  Y  are  known  to 
this  function  through  a  common  block  named  DATA4.  So  only  NPNTS1, 
NPNTS2,  NPNTS3,  IP0WR1  (-  k-1),  IP0WR2  (-  L-l),  IP0WR3  (-  m-1)  and 
IY  (-0)  have  to  be  passed  to  the  function  and  the  value  of  is 

returned  through  the  function  name  S4. 

** 

SUBROUTINE  MAT42  presets  the  subaatrices  2.  class  A.  in  the 

^*1 

following  way 

i)  Preset  edge  section  elements  (using  REAL  FUNCTION  S4) 


./ 


,oa  ,  ,o... 


e^ak;liN 

G 

a*;Lj  V*/ 
%l;  Ht2 

4 


~~  A*}h 
tfjU  -me 


$ 


ii)  Copy  defined  elements  diagonally 


•  00 
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*>S? 


»,v 


J,V 


^k;/;  V,/ 

flW;V  ^kyijU  a*;/;/.J  *■•  aU/./» 

^LV  %;/.2  ^.US  •••  S./.V* 

°V^g> 


a";tiU 

akU,U 

\/;<7 

°*i 1:22 

**&•« 

°*iU.I 

auv 

•  •  • 

•  •  • 

#  •  • 

Vv 

v« 

^UfU 

‘V*' 

•  •  • 

ak;l;A2 

•  •• 

%Ut 

•  •  • 

aH)K> 

®u*u 

k,h*JK 

ak;l:  /,M 

aTeiLii-M 


*.  1;  ',M  I 


V;*," 


hhl* 


Vi/  KM  1 


72 


k 


A 


* 

k 


A  u 

A** 

<r 

s« 

*  ** 

a*' 

t'kU 

tiu 

**• 

A** 

> 

/  * 

/* 

a#  • 

4** 

'"A  >*./» 

Kl 

KZ 

Ak;u 

A** 

•  *  *  /,#• 

*Vu 

*  e-P 

"V-Zj 

A** 

... 

A4* 

•••  “ksJjf 

••  •  —  •  •  •  .  #•  «  a 

/S  S  eS 

•  •  •  »  •  •  • 

Akntf 

.  ** 

A** 

\>.7 

A** 

>** 

*k,2j 

<1 

A** 

A** 

••  • 

A** 

*«0 1' 

K.u 

...  Ak  ;j,n 

•  •  0 

A** 

•  •  • 

c. 

•  •  • 

A*1 

e«w  • • • 

•  •  •  A  k;M,n 

This  subroutine,  that  naads  NPNTS1 , NPNTS 2 , NPNTS3 ,M,  N , IPOWR1 ("L-l ) 
as  Input  paraaatars,  returns  the  subroutine  1st  class  a£  (SUBM1(16,16)) 
SUBROUTINE  MAT4,  which  presets  the  entire  systea  matrix  A  in  the 


following  way. 


to 


J 
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4.4.2.  Right  Hand  Side  Vector  B: 


The  right  hand  side  subvectors  1.  class  generally  can  be 


written  as 


mvC 

2 2  EY„-x£xjrt 

•  e  e 

2  2  2  VC*£ 

22  2t-<-x^ 

2  2  SY„{Xi^')f2»cX?lli 


EEZY^'^-Xj*' 


2  2  2X,« 

2E2x<</*C-x?* 

2  2  2  Y^'Xli'-xC-xC 


k*  /. ... .  i 
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The  subveccors  1.  class  (k 
up  into  subvectors  2.  class: 

B** 


1,...,L)  themselves  will  be  divided 


v 


h,L 


where  subvectors  2.  class 


** 


p** 


B  ^  can  be  written  as 

ESeVv<!-X2|2 

x3 


H-l 


ESEXsKl-nZ-x* 

k-  !j  ...,  L  j  L*  /y  ...  ,  M  ;  «'»*  /,  .. 


A/ 


Finally,  the  elements  of  the  subvectors  2.  class  ^ 

h  -  ZllY^Xl"-XZ"-X3r' 

K)  L,  m 


can  be  written  as 


1  «l  wt 
Lx/  J  .  L  *  !  M  • 

"  /•••/“*  x  J  •••  J  / 


’r* 


m  *  fj  ...  y  N 


To  calculate  b^^,  REAL  FUNCTION  S4  is  used,  again.  Data  XI,  X2,  X3 
and  T  are  available  through  COMMON  block  DATA4 .  NFNTS1,  NPNTS2,  NPNTS3, 
IP0WR1  («  k-l),  IP0WR2  (-  L-l)  EP0WR3  (-  ®-l)  and  IT  (-1)  have  to  be 
passed  to  the  function  and  the  value  of  bk.^.a  is  returned  through  the 
function  name  S4.  Using  FORTRAN  programming  language,  the  allocation 


of  b.  can  be  performed  by  three  stacked  DO- loops. 

f® 
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DO  07  U-1,L,1 

IPOWR1-L1-1 

DO  08  M1-1,M,1 

IPOWR2-M1-1 

DO  08  Nl-l.N.l 

IPOWR3-N1-1 

1*1+1 

08  B ( I ) -S A ( NPNTS 1 , NPNTS  2 , NPNTS 3 , IPOWR1 , IP0WR2 ,  IPCWR3 , 1 ) 


4.5.  Software: 


The  software  to  compute  the  coefficients  for  a  two  dimensional 
approximation  is  described  in  APPENDIX  A  and  is  Implemented  in  the  TPL 
HP-21MX  computer  system.  To  work  correctly  with  these  program  modules, 
the  user  has  to  conform  to  the  following  conventions. 

i)  Provide  Che  data  in  four  arrays 

(Type:  REAL)  of  5*5*5*  elements 

through  a  COMMON  block,  named  DATA4. 

COMMON  /  DATA4  /  X1.X2.X3,? 

REAL  n(5,5,5),X2(5,5,5),X3(5,5,5),Y(5,5,5) 

ii)  Dimension  an  array  (Type:  REAL)  of 

5*5*5  elements  to  contain  the  coefficients, 

REAL  COEP  (5,5,5) 

ill)  Define  the  parameters  NPNTS1,  NPNTS2, 

NPNTS3,L,M.N  and  IPR2NT 

NPNTS1...  #  of  XI  variations 
L  i  NPNTS1  i  5 

NPNTS2...  t  of  X2  variations 
M  i  HP NTS 2  i  5 

NPNTS3 . . .  #  of  X3  variations 

N  <  NPNTS3  i  5 

L  ...  (desired  order  of  approximation 
polynomlnal  v.r.t.  XI)  +•  1 
1  <  L  i  4 

M  ...  (desired  order  of  approximation 
polynomlnal  v.r.t.  X2)  +  1 
1  (  N  (  4 

N  ...  (desired  order  of  approximation 
polynomlnal  v.r.t.  X3)  +  1 
10(4 


IPRINT  ...  Controls  quantity  of  print  out 

2  ...  Print  system  matrix  and  right 
hand  side  vector  before  and 
after  Gauss  Jordan  Elimination. 

1  ...  Print  equation  system  after 
Gauss  Jordan  Elimination 
<0...  No  print  out 
>0...  Print  equation  (4.1.)  with 
the  actual  parameters. 

iv)  If  a  user  program  uses  subroutine  MAT4,  the  software 
modules  have  to  be  loaded  using  the  procedure  in 
section  1.4 .iv. 

If  all  these  requirements  are  met,  the  correct  call  for  the  subroutine  is: 

CALL  MAT 4  (NPNTS1 , NPNTS  2 , NPNTS  3 , L , M , N , COEF , IPRINT ) 


Upon  completed  execution  of  this  approximation  routine  array  COEP 
contains  the  coefficients.  Externals  used  by  MAT4:  AB4,  DATA4 ,  MAT41, 
MAT42,  S4.  Under  no  circumstances  may  the  user  use  any  of  these  names  for 
a  modules  of  his  own  user  program. 
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4.6.1  PTN4  Compiler  Hating  of  aamplo  user  prograa 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

ms 

0013 

0014 

0015 

0016 

5317 
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FTN4.L 


C 

C 

C 

C 

C 

c 

c 


POOCH AH  FITR  <3,99) 

APPROXIMATION  POOCH AH  TO  FIT  A  SURFACE  THROUGH  DATA  POINTS  Y , 
WHICH  DEPEND  ON  THREE  INDEPENDENT  PARAMETERS  XI,  X2  AND  X3. 

AUTHOR:  HANS  ZEBNER 
DATE  :  AUGUST  26,  1980 

ft,  Han*  Z«bn«r:  4D-Appr  oxidation*,  Y«T(X1 ,X2,X3) . 

COMMON  /  DATA4  /  X1,X2,X3,Y 
COMMON  /  AB4  /  A,i  ’ 

REAL  Xl<5,5,5) ,X2(  5 ,S,S) , X3< S , 5 , 5 ) ,Y(5,5,5) 

REAL  A<64,64) , B(64> 

, IFILE<3> ,N0CR(2) ,ICLR<3) 

DATA  NOCR  /000033B , 040433B/ ,  ICLR  /015S24B, 01SS1SD, 006S37B/ 

2A2 


101  FORMAT  <■  ENTER  NPNTS1  NPNTS2  NPNTS3 
*) 

<•  ENTER  L  M  N  -2A2) 

(•  ENTER  DATA  FILE  NAME  *2A2> 

<3A2) 

<*  ENTER  IPRINT  »2A2) 

FORMAT  <*  OUTPUT  RAW  DATA  YES  OR  NO  a2A2> 
FORMAT  </•  XI •) 

FORMAT  </lX,I2" . 


102  FORMAT 

103  FORMAT 

104  FORMAT 

131  FORMAT 
13: - 

106  FORMAT  (>1X 

107  FORMAT  “ 

108  FORMAT 

109  FORMAT  .. 

110  FORMAT  </*  X 


DATA  SET-/) 
14F10.3) 


<3X,i2I10> 

< IX* 12, 14F 
(/•  X2*) 

...  ...  </•  X3-) 

lit  FORMAT  </"  Y") 

til!  \wmv 

READ  RAH  DATA  FROM  FILE 

01  WRITE  (1,101)  NOCR 

READ  <1,  *>  NPNT31 , NPNTS2 , NPNTS3 

WRITE  (1,149)  ICLR 

NMAX-S 

» a  a  w 

IF  (NPNTS3.GT.NMAX)  CO  TO  01 

'  ,102)  I  "■ 


02  WRITE  (1, 

READ  (i;  ft) 
WRITE  U^IA*) 


NOCR 


IF  (ttSMftN.GT.64)  GO  TO  02 

. '  WSp 


.  (LftMftN _  _ 

03  WilTE  (1,103)  N00R 


MW 


nil 

0060 


s  m 

h# 

m 

0071 

0072 

0073 

to75 

0076 


WIT* 

31  fill*  IPRINT 

Bill  MW 

IL-2*NPNT81*NPNTS2ftNPNTS3 

CALL  READ?  ( IDCB .IERR, XI , IL, LEN, 1 ) 

if  (IERR.lt.o)  stop  0601 

^VI||g0^^s?clt^"o6gi'IL'LEN’S4, 

S^i|ISD[T<I?c!tiP?igi-I,-'LEN'107 

CALL  READ?  ( IDCB , IERR . Y . IL, 

IF  (IERR.LT.O)  STOP  000* 

CALL  CLOSE  (IDCB. IERR) 

IF  (lERR.LT.n)  STOP  0005 

tfilTOl*?'"6’  co  TQ 17 


107) 
LEN, 160) 
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PACE  0002  FITR  3 i 47  PM  FRI.,  26  SEP.,  I960 


Results 

printed  output 


XI 

1. 

1 

I 

4 

5 


DATA  SET 

t 

.  100 
.  100 
.100 
.  too 
.100 


2 

.100 

.100 

.100 

.100 

.100 


2.  DATA  SET 


3 

.100 
.  100 
.  100 
I  100 
.100 


4 

.100 
.  100 
.100 
.  100 
.100 


s 

.100 
.100 
.100 
.i  eg 
.100 


1 

1 

4 

s 

3.  DATA  SET 


i 

4 

5 


4. 


.  20? 

.m 

'.loo 


.20$ 

.loo 


4 

.200 


>00 

Soo 


.20? 

'.ill 

.loo 


1  .40? 

?  .400 

3  .400 

i  -.m 

S.  DATA  SET 

.40$ 

.400 

,400 

Ml 

.40$ 

All 

'.ill 

.40$ 

.400 

.400 

All 

.40? 

Ml 

:ttl 

1  .SOO 

£  .soo 

3  .soo 

4  .500 

s  .soo 

X2 

2 

.SOO 

.soo 

.soo 

.soo 

.soo 

3 

.soo 

.'ioo 

.500 

4 

.500 

:!it 

.500 

.500 

.50? 

:lo! 

.500 

.SOO 

1.  DATA  SET 

1  -20.00? 

2  -10.000 

I  iiilli 

2.  DATA  SET 

-20.00$ 

-10.000 

liiiil 

-20.00$ 

-10.000 

iiiSSS 

-20.00$ 

-10.000 

0.00§ 

mn 

-20.00? 

-10.000 

0.000 

mn 

1  =**:Ki 

1  .t.-tii 

S  20.000 

3.  DATA  SET 

2 

20.000 

'IMS! 

20.000 

3 

-20.000 

'ISiSSS 

20.000 

4 

-20.000 

'IMS! 

20.000 

-20.00? 

•ISiSSS 

20.000 

i 

1  48: SIS 

4.  DATA  SET 

-20.00$ 

11:111 

18:888 

-20.000 

■T.m 

mn 

-20.00?. 

"l?:m 

mn 

1 

2 

3 

4 

s 
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-20.000 
-10 . 000 
0.000 
10.000 
20.000 

data  set 

-so.oolk 
-10.000 
0. 008 
10.000 
20.000 


1.  DATA  SET 

1  -20 . 00^ 
2  -20 .008 

3  -20.800 

4  -go. ooo 
s  -20.000 

2.  DATA  SET 

1  -20 . 00^ 
2  -20.000 
3  -go. 000 

i  --inn 

3.  DATA  SET 

i  :|S.SIS 
I  -IS: III 

5  -20.000 

4.  DATA  SET 

1  -20.0ol 

2  -20.000 

3  -20.000 

4  -20 .600 

5  -20.000 

5.  DATA  SET 

i  -IS: So. 

3  -20.000 

4  -20.008 

5  -20.000 


-20 . 000 
-10  .000 
0.000 
10.000 
20.000 


-20.000 
-10 . 000 
0  .  000 
10.000 
20.000 


-20.000 
-10 . 000 
0.000 
10.000 
20.000 


-20.000 
-10.000 
0  .  000 

UMl 


-20.  oof 
-10.000 
0.09a 
to.ooe 
20.000 

-20. oof 
-10.000 
@.0  00 

mu 

-20 . 00$ 
-10.000 
0.000 

mn 

• 

-20. oof 
-10.000 
0.000 

mn 

-iO.Oof 

-10.000 

-10.000 

-10.006 

-10.000 

o.oof 

0.000 

0.000 

0  .  0  00 
0.000 

10.000 

10.000 

10.000 

10.000 

10.000 

20.000 

20.000 

20.000 

20.000 

20.000 

-lo.ooi 

-10.000 

-10.000 

=it:!SI 

o.oof 
0.000 
§ .  000 

S:IIS 

10. oof 
10.000 
10.000 

mn 

20. oof 
20.000 
20.000 

inn 

-10.000 

--mn 

-10.000 

-10.000 

o.oof 

0.000 

0.000 

0.000 

0.000 

10.000 

mn 

10.000 

10.000 

20.000 

inn 

20.00c 

20.000 

-10. oof 
-10.000 
-10.000 
-10.000 
-10.000 

o.oof 

mi 

0.000 

0.000 

10.000 

mn 

10.000 

10.000 

20. oof 

mn 

inn 

-10. oof 
-10.000 
-10.000 
-10.000 
-10.000 

o.oof 

0.006 

0.000 

0.000 

0.000 

10.000 

10.000 

10.000 

10.000 

10.000 

20. oof 
20.000 
20.000 

inn 

1.  DATA  SET 


2.  DATA  SET 


sof 

.700 

.  4  of 

Iff 

:!ff 

m 

:?f°0 

Ml 

4 

5 


1 

5 

4 

5 


i 

4 

5 


.463 

DATA  SET 

.763 

.863 

.763 

.463 

1 

.650 

•Ml 

.700 

.650 

DATA  SET 

2 

.  950 
1.000 
1.050 
1.008 
.950 

3 

1.050 

1.100 

1 .  ISO 

i:^S 

.950 

1:828 

«6sl 

:Jlt 

'.600 

1 

.963 

1.013 

1.063 

1.013 

.963 

data  SET 

2 

1.263 

m 

1.263 

3 

1.363 

1.413 

1 . 463 
1.413 
1.363 

4 

1.263 

m 

1.263 

5 

.963 

1:113 

Ml 

i.ui 

1.500 

1.450 

1.400 

Mil 

1.7S0 

1.700 

(«l 

1.850 

1.800 

i  ai 

i:?S 

(ill 

im 

1  i  I  :• 

K  •  i 

i  -a 

1  ■  i  .4S8e"o1  106£+9§ 

5  =  5  -:135iIS3  : 3111:85 

k  »  i  a 

I  ■  3 

/  •  I  -.994E+0S  . 230E+05 

} :  I 


■.7S7E-.3 

: 1251:81 

3 


> 16SE+04 


K  -  1 


-.41 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Compared  Co  che  probe  calibration  approach  described  in  Ref  3. 

Che  method  presented  here  provides  a  much  simpler  way  to  both  calculate 
che  probe  calibration  surfaces  from  measured  data  and  to  apply  the 
calibration  to  on-line  data  reduction.  The  iteration  required  in  the 
method  of  Ref  3  is  completely  omitted.  It  is  noted  however,  that  the 
Gauss  Jordan  Elimination  lacks  some  sophistication.  The  system  matrix 
for  example  undergoes  the  elimination  without  prior  conditioning. 
Particularly  in  the  case  of  the  3-D  and  the  4-D  approximations,  numerical 
round  off  errors  have  an  influence  on  the  accurraey  of  the  coefficients. 
Consequently  it  is  recommended,  that  the  Gauss  Jordan  Elimination 
routine  be  revised  to  use  double  precision  constants  and  a  routine  that 
conditions  the  system  matrix  be  added.  Time  constraints  made  these 
steps  impossible  for  the  author. 

Following  these  refinements,  a  calibration  procedure ,  similar  to 
McGuire's  (ref  3)  should  be  formalised.  Since  che  4-D-  approximation 
requires  large  arrays  and  therefore  extensive  CP  -  memory,  extended  memory 
access  (EMA)  is  necessary  if  the  H?  21MX  computer  is  used  to  perform 


the  calculations. 
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Some  Useful  Matrix  Cooventiona  and  < 
Submatrix  Notation 


az, 

<*4 

*A, 

*44 

APPENDIX  B:  Software  Description:  Flow  Charts 


The  following  are  given: 


Flow  Chart  Page 

MAT2  93 

52  97 

MAT 3  98 

MAT31  106 

53  109 

MAT4  111 

MAI41  119 

MAT42  124 

54  127 
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Flow  c-l^ort 


^A0-A110  673  BOH  CORP  MONTEREY  CA  ...  . 

rsrsrs  jSncT^  p*o"*,,  ro*  ?*  »  BaT/m?(u> 


UNCLASSIFIED 


00  04  fC  / 


\ 

yeoQe  eeebe*  now 
s&»s  pmaeU 


/W  QZ 


Define  location 

Define  locator) 

of  element 

of  element 

Tel 

r-L 

JmK 

Gonyou**  tSpuotA  'on 

2  +•  ^"y 

Furicnorr  SZ 
and  write  value 
Info  A  (r,  JJ 


Proceed  -fa  next 
edje  element . 
Increase a outer. 
IPotue/  •IFotdfi*/ 


“ere  aU  \ 

lSw?7/ 


94 


<shari  MAT3 


■Dtfir*  banh'ato 
of ssubmcrtrix 
l  s/or/ 

/  s/op  c/m+n/s 
TSTAfiT.l 
JS70P-  M 
JSTART*  OC-0*ft+! 

J srw*-  K*M 


Define  loco /ion 
of  suitmntri'r 
/  doss  bysiarl 

$  s^odn*n*/s 
ISTAR  T*  f K-U*tf+  / 
J&roA  {CK-*  I) 
JSTAKT*  a_-l)*n+l 
JSrop-  L  *M 


>o* 


100 


4 


Hoce*e(  /a 

Increase  po**r 

J&ut/r/  *&x*tt+( 


/4r*  aO 

44jm  ■»***  * 


Ccf>y  cf*f<ned 
dtrmw'h 
«/»ay*#W/p 
U-L-t 


bo  or  D’l.U.l 


a>drjt*a,L,/ 


simri  i 

rleft  (jbiNn/f  fir 

**  «vy 

JSTAfT.  t 

ZW»*  -Z7*^f 

Jsnor  (?i-/)*m*i 
Osror  *  7/«^ 


DO  Oi  I'KTAKT.JSTCF,  I 

Do  <*  3'&rMr,3rra°,  i 


104 


0 


i 


equation 

(T.7. 

e/iry»t*a.£>o)n 

/e\ 

_ i 

L 

r 

H'/rit  the 
coettc/err/s 

. — 1 

1 

RETURN 


ow 


c Inarl  ^AT^f 


(return  ) 


108 


/7ow  chart  <S<3 


/ 


fil 


~/ow  chari  MA  T4 


J.12 


DO  04-  t-TSTA*Tjrro*,) 


M 


2  *I/+  ' 

J!'  0 


Do  o*  J- jsrrxrr,  xro*.  / 


J/'J/*/ 


VI 


Copy  -tuhrminx 
SU8MI  at  dmr,  vW 
wkrouiiHt 
MATH  mb  aomeci 
fib  C*  ir>  A 
AfIJ)'OJtHtaS'3) 


m  Am  tv*r 
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Recced  <o  noxi 
'  i-Amahn* . 

,jV»  r-'^se  po %*cr. 

Zfajent 


Mlf  «*✓  \ 

•tUe  suh^^c 
_  preset? 


Coftj  defined 
SuemaSncer 

ckajono/(y 
U-L -/ 


’  J7-  /,  *./,  / 


Jo  o7  4.,  / 


Oefee  stmr4  em4 
Shaft  s/en*»**  far 

Hm  copy 
JSTMr* 
lsn^»n*MN 
zsmrr-  cjr-/)*^mi 

J57&n  CM/W 
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CO  fotar 
Compietod 


PtiM  equ  oJien 
Sy*hm.m  4. 1. ; 
Sfsitm  matrix  A 
KmM  itand  c id*  1*0 4ri 

ZP*//vr*7 


Gauss  Jordan 
oknnafyort  wdh 
matrix  A  A 
vector  H 
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1 


tyuaAot* 
oflebG.  7 
o/ithnnot/iotn 

/ a 

\ 

u, 

r 

/W  A 

G>*ffic\ 

w  1 

en/jr  1 

_ s 

1 

{^return  ") 

"W  chad  MAT41 


(ENTRY  ) 


io  pr&ei 
edge  sector, 
suAnnaAttcet 
Z.  c/ass 
TP0W*2  *  o 


Do  o4-  K*I.MI,i 


D*Ano.  /ooadon  of 
SoAmaJnr  Z.  dan 
Ay  dot*  4  tf/eye 
derntnfa 

TSTWl 

lsmr*N 

3srAgr*0(-t)*n+t 


'Afioe  /oca  Aon  of" 
Sobmaihi*  2.  don 
Ay  Siar6  4  s(y» 
demonts 
SSTWr*  0<-M)+t*+l 
P*ve  *<Hr*i)-n  )*n 

OSTOAt 


04 

OCO 


i 


Preset  edjr  socAon 
SuAmafn*  2.  c/atr 
/$u 9tm )  uimj 
****urm  Hath. 
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© 


00  aSI-Tnjmrjsroe.i 


X1-B+/ 

?/• 0 


jo  aJ  J-Jstakt,  zsror,  / 


_ 9 

Cafij  SuitmQtnr 
?.  cAuv  SV&T2 
Of  derive*  fiwn« 
fl/hroufane  H4T*2 
into  correct  pimct 
in  aukmmfrt* 

I.  door  SutM/. 
Sttt»rr,3hajiHt/e&) 


Orfine  siari  omd 
step  afe 
for  Ma  copy 
Ttr*rr.  (zh/)*N+i 

zsjoo*  zt*n 

JStMT-  (Oh  !)•**> 
J SToh»J/*H 


Oogti>lsmr,jsroP,  / 


Do  os  jrjsnor,Dsno,  / 


1 


<2}»v 

•Mwa*,  jw*  jww  rr.j; 


Flow  chari  MAF4Z 
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APPEND1"!  C :  Software  Description;  Listings 


The  following  are  given: 


Listing 

Page 

Two-Dimensional  Approximation 

132 

BLOCK  DATA  AB2 

132 

BLOCK  DATA  DTA2 

133 

PROGRAM  DEM02 

134 

SUBROUTINE  MAT 2 

137 

REAL  FUNCTION  S2 

140 

INTEGER  FUNCTION  IEL22 

141 

SUBROUTINE  IEL2 

142 

REAL  FUNCTION  F2 

145 

Three-Dimensional  Approximation 

146 

BLOCK  DATA  AB3 

146 

BLOCK  DATA  DATA3 

147 

PROGRAM  DEM03 

148 

SUBROUTINE  MAT3 

151 

SUBROUTINE  MAX31 

155 

REAL  FUNCTION  S3 

156 

INTEGER  FUNCTION  IEL33 

157 

SUBROUTINE  IEL3 

158 

REAL  FUNCTION  F3 

161 

Four-Dimensional  Approximation 

162 

BLOCK  DATA  ABA 

162 

BLOCK  DATA  DATA* 

163 

130 


%wV*»' 


APPENDIX  C  (Cont'd) 


Four-Dlaenaional  Approximation  (Cont'd) 


PROGRAM  DEM04 

164 

SUBROUTINE  MATA 

168 

SUBROUTINE  MAT 41 

172 

SUBROUTINE  MAT42 

173 

REAL  FUNCTION  S4 

174 

INTEGER  FUNCTION  IEL44 

175 

SUBROUTINE  IEL4 

176 

REAL  FUNCTION  F4 

179 

131 


far1** 


MNI4N  Otlflt 


Man# 


Ml  «••#*#• 


I 

\ 


I 

r 

i 


l 

s 

ti 

n 

ti 

n 

M 

ss 


*1 

X 

99 


•9 

ft 

T» 


•• 


•W/ll/2*a  l*a«Ca*l 


^•oiiu  omoj  <i»«jT»jur*jf«rA*ii«rAr<o*rAM«*(n»uri 

MaiOttittMiMMiiaitmiiiMtaMiMitiutMuititaMMiMaMiai 

:  B'USk'SJKB'ttiR.rSJM  «&.«,?:*•*  u‘  e9,*,c,  : 
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